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ABSTRACT 

Radio and optical images of the M87 jet show bright filaments, twisted into an 
apparent double helix, extending from HST-1 to knot A. Proper motions within the jet 
suggest a decelerating jet flow passing through a slower, accelerating wave pattern. We 
use these observations to develop a mass and energy flux conserving model describing 
the jet flow and conditions along the jet. Our model requires the jet to be an internally 
hot, but subrelativistic plasma, from HST-1 to knot A. Subsequently we assume that 
the jet is in pressure balance with an external cocoon and we determine the cocoon 
conditions required if the twisted filaments are the result of the Kelvin-Helmholtz (KH) 
unstable elliptical mode. We find that the cocoon must be cooler than the jet at HST-1 
but must be about as hot as the jet at knot A. Under these conditions we find that 
the observed filament wavelength is near the elliptical mode maximum growth rate and 
growth is rapid enough for the filaments to develop and saturate well before HST-1. We 
generate a pseudo-synchrotron image of a model jet carrying a combination of normal 
modes of the KH instability. The pseudo-synchrotron image of the jet reveals: (1) that a 
slow decline in the model jet's surface brightness is still about five times faster than the 
real jet; (2) that KH produced dual helically twisted filaments can appear qualitatively 
similar to those on the real jet if any helical perturbation to the jet is very small or 
nonexistent inside knot A; (3) that the knots in the real jet cannot be associated with 
the twisted filamentary features and are unlikely to be the result of a KH instability. 
The existence of the knots in the real jet, the limb brightening of the real jet in the radio, 
and the slower decline of the surface brightness of the real jet indicate that additional 
processes — such as unsteady jet flow and internal particle acceleration — are occurring 
within the jet. Disruption of the real jet beyond knot A by KH instability is consistent 
with the jet and cocoon conditions we find at knot A. 

Subject headings: galaxies: individual (M87) — galaxies: jets — galaxies: active - 
radio continuum: galaxies — hydrodynamics — relativity 

1. Introduction 

M87 (Virgo A, NGC4486, 3C274) is a giant elliptical galaxy near the center of the Virgo 
Cluster. The well-known, kpc-scale jet in this galaxy is a prominent source of radio, optical, near- 
ultraviolet (NUV) and X-ray emission. There are remarkable similarities between the radio, optical 
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and NUV emission on scales ~ 0.1 — 1 kpc. The bright knots, the filamentary emission between the 
knots, and the bends and twists in the jet can easily be identified in the optical and NUV as well 
as the radio. While X-ray images from Chandra are not of comparable resolution to the radio or 
optical, the same bright knots and overall jet structure can easily be seen in the X-rays. Relativistic 
proper motions are seen within the knots, both in the radio and the optical. 

A variety of models have been proposed fo r the p hysical state of the kpc-scale jet, but none 
has emerged as definitive. iBicknell <fe Begelman (ll996l ) proposed a hydr odynamic model, in which 
the knots are shocks caused by the Kelvin-Helmholtz (KH) instability. iHeinz fc BegelmaiJ (119971 ) 
extended the hydrodynamic model, arguing that intensity changes and spectral evolution along 
the jet can be explained by shock compression of a weakly magnetiz ed plasma. Other author s 
have put forth models in wh ich the magnetic field plays a key role. iFraix-Burnet et al.l (jl989l ) 
and also IWalker et all ( 19941 ) proposed that the jet is a force- free MHD configuration (following 



Konigl k, Choudhuril (|1985l )). I n this scenar i o the knots correspond to "mag netic islands" in the 
force- free state. More recently, iGracia et al. ( 20051 ). also iGracia et al.l (|2009h . presented an MHD 



model of the jet's launching from the central black hole; they identified HST-1 as a recollimation 
shock. 

In this paper we return to hydrodynamic models and consider the pair of twisted emission 
filaments reported by Lobanov, Hardee & Eilek (2003; hereafter LHE). We propose that these fila- 
ments are created by the KH instability in a weakly magnetized, predominantly hydrodynamic flow. 
Several authors have used the properties of KH modes to estimate the physical conditions in jets 
and their surroundings (e.g., Perucho et al. 2004). The method has been applie d to knots in the th e 
M87 jet (IBicknell fc Begelman|[l99^ ): to twist ed structures in the 3C 120 jet (lHardee et al.ll2005l ): 
to twisted emission threads in the 3C 273 jet (jPerucho et al J 120061 ) . to super lumina l motions and 
accel erations of components along curved trajectori es in the 3C345 iet (lHardedll987 ; Steffen et al.l 
19951 ). and to tr ans verse jet structure in 0836+710 (jPerucho fc Lobanovll20o"7f r 

We bas e our study of the M87 jet on a linear analysis of the KH instability, e.g., lHarded 
(120001 . 120071 ). supplemented by our experience with numerical simulations of its development and 
saturation. We chose this approach because our goal is to constrain the physical conditions in 
the M 87 jet and its surroundings, and it would be prohibitively expensive to run the large suite 
of simulations that would be necessary for a good exploration of parameter space. For example, 
computational constraints make it difficult to simulate an expanding, three-dimensional, relativistic 
jet, close to the line of sight, over the distance from HST-1 to knot A, with sufficient spatial 
resolution and temporal storage. Such a simulation woul d need to have enou gh computational 



zones to satisfy projection and numerical viscosity issues ( Perucho et al. 20041 ) . The simulation 



would also have to span and store enough time i nformation to deal with the light-tra vel time effects 
that lead to the observed superluminal motions ( Gomez et al. 1997 ; Alov et al.ll2003l ) and that make 



the ap parent structure different from the intrinsic structure. However, more recently iMimica et al 
(120091 ) have developed a Lagrangian algorithm that accomplishes this task in a better and more 



physically correct way. 

We proceed and organize the paper as follows. In Section [2] we review the structure of the 
observed jet, particularly the filaments identified by LHE, and the proper motions of the bright 
knots. We use the proper motions to estimate the flow speed of the jet plasma and the wave speed 
of the KH pattern. The data suggest that the jet decelerates between HST-1 and knot A; this is 
a key point of our models. In Section [3] we develop an energy-conserving jet model which matches 
the observed flow deceleration and determines the thermal state of the jet plasma. We also develop 
models in which a small part of the energy flux is lost to the cocoon. In both models we assume the 
jet is in pressure balance with an external "cocoon" , as required by our KH analysis. In Section H] 
we present a KH instability analysis, and determine what the thermal state of the cocoon must be 
if the elliptical mode of the KH instability causes the observed filaments. We do this by matching 
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solutions of the wave dispersion relation to observations of the pattern speed and wavelength along 
the jet. We also verify that the KH instability growth rate under these conditions is sufficiently 
rapid for the instability to saturate. In Section [5] we explore the likely internal state of the jet and 
cocoon plasmas. In Section [6] we present a pseudo-synchrotron image of the jet, to compare our 
most likely model to the observed intensity structure of the real jet. Finally, in Section [7] we briefly 
summarize our results and in Section [8] we discuss a few implications of our models. 

2. Jet Structure and Environment 

In this section we review the basic structure and immediate environment of the M87 jet. We 
also review low-order wave modes created by the KH instability and compare them to the structure 
of the real jet. Because our KH analysis requires a model of the jet flow speed as well as the speed 
of the wave pattern within the jet, we use observed proper motions to determine the evolution of 
both speeds along the jet. 

The overall structure of the M 87 jet as seen in the radio is shown in Figure [1] (from Owen, 
Hardee & Cornwell (1989); hereafter OHC). The brightest emission knots in the outer part of the 
jet are labelled A, B and C. The inner jet — from the core to knot A — is expanding uniformly, 
limb brightened with sharp edges. Bright knots and fainter filaments can be seen within the flow. 
Past knot A the jet brightens, stops expanding, and maintains a more uniform radius. Knots B 
and C contain tightly twisted filaments as well as some bright emission patches. Past knot A the 
flow develops significant sideways motion, and past knot C it breaks up dramatically. We take the 




core 



Inner halo 



Fig. 1. — Radio image of the jet and western part of the inner radio halo in M87. The three bright knots A, B, C 
are labelled, as is the western part of the inner radio halo. The straight, well-collimated inner jet, which we model, 
extends from the core to knot A. The jet begins to destabilize past knot A, and disrupts fully past knot C, where the 
plasma flow bends and feeds the inner halo. VLA image taken at 15 GHz, with resolution 0.15 arcsec, from OHC. 



distance to M87 to be 16.7 Mpc (|Mei et al.ll2007l ). thus 1" ~ 80 pc. As we discuss below, the jet 
probably lies between 15° and 25° from the line of sight, so its intrinsic length is ~ 2.4—3.9 times 
longer than its apparent (projected) length. 

The plasma content of the jet remains unknown. Because the jet is a synchrotron source, we 
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know it must contain a relativistic, magnetized lepton plasma. Reynolds et al. ( 19961 ) argued for a 
pure lepton jet, because a cold, electron-ion jet would violate constraints on the opacity and surface 
brightness of the self-absorbed, pc-scale core. However, mor e recent evidence on the je t viewing 
angle (as in Section I2.3j) and the internal energy of the jet ( Bicknell Sz Begelmanl 19961 ) and our 
results in Se ction 13.31 expand the p ossible parameter space to allow electron- ion jets within the 
arguments of lRevnolds et all ( 19961 ). 

The local environment of the jet, on projected scales £ 2 kpc, is well studied in radio and 
X-rays. From the radio, we know that the jet feeds, and appears to be surrounded by, inner radio 
lobes which extend ~ 3 kpc from the galactic core. The western part of the inner lobes can be seen 
in Figure [TJ The inner radio lobes contain a relativistic, magnetized lepton plasma. The minimum 
pressure allowed by the synchrotron power of the inner lobes is below that of the local interstellar 
medium (ISM) th roughout most of the lobes, but a few bright filaments are at higher pressure 
( Hines et al. 19891 ). Such bright filaments are probably transient, suggesting jet-driven turbulence 
in the radio lobes. 

From the X-rays, we know that the jet sits within the central ISM of M87. That ISM is an X- 
ray loud, ion-electron plasma; a verag ed over the inner f ew kpc it has T ~ 10 7 K and n <y 0.1 cm 



e.g-.lNulsen fc Bohringerl ( 19951 ). also lOwen et al.l ( 2000l ). The inner ISM is magnetized; lOwen et~al 



( 1990l ) detected Faraday rotation from the ISM in front of the jet and inner radio lobes. The inner 



ISM is disordered a nd turbulent; filaments, lo ops and bubbles are apparent in the inner few kpc of 
the Chandra image ( For man et al. 2005 . 20071 ). The correlation of features in the X-ray image with 
those in the radio halo shows t hat the rad i o-loud plasma from the jet is interacting, and probably 
mixing, with the thermal ISM. Keel et al. (j 19961 ) found that emission-line clouds in the inner few 
kpc of M 87 show disordered motion at about half the sound speed of the local ISM; these clouds 
probably share the velocity field of the ISM. 

2.1. Twisted Filaments in the M 87 Jet 

Figure [2] shows the inner jet in M 87. The prominent emission knots are labelled, from the 
core outwards, as HST-1 (formerly called knot G), D, E, F, I and A. HST-1 is ~ 0.9" (projected 
distance ~ 70 pc) from the core; knot A is ~ 12" (projected distance ~ 1 kpc) from the core. 
The jet expands with a constant apparent opening angle, ~ 6°, from well inside HST-1 to knot A. 
Between HST-1 and knot A the emission profile is edge-brightened, relatively sharp in the radio, 
and contains the twisted filaments reported by LHE. The projected magnetic field ve ctors lie mor e 
or less along the edge of the jet flow (OHC), which may be indicative of a shear layer ( Laingj[l98ll ). 
However, there is no i ndication of mass ent rainment into a turbulent surface layer as suggested for 
3C31 on large scales ( Laing &: Bridle! 120021 ). 

The jet structure seen in the optical and NUV is quite similar to that seen in radio, but there 
are important differences in detail. The optical/NUV emission is more concentrated in the knots 
and towards the Jet axi s, and the outer edges of the jet are less well defined than in the radio 
e.g.. ISparks et al.1 (Il996l ) or lMadrid et al.l ( 20071). The o ptical/NUV knots also differ somewhat in 
polarization from the radio knots ( Perlman et al. 19991 ) . The major bright radio or op tical knots 
can also be identified in the lower-resolution X-ray image, e.g., iPerlman fe Wilsonl (|2005h . but they 
can differ in position or structure from their radio or optical counterparts. 

In this paper we do not focus on the knots, but rather focus on the twisted filaments which 
LHE detected in the inner jet of M87. LHE extracted slices transverse to the overall jet axis, 
spaced by the pixel size, in radio (VLA) and optical (HST) images. They fit each slice with double 
Gaussian profiles, allowing the amplitude, position and width of the two Gaussians to vary at each 
slice. Both VLA and HST images revealed a consistent pattern: the jet contains two intertwined 
emission filaments, which can be traced from HST-1 out to knot A. Unlike the situation with the 



Fig. 2. — The inner part of the M87 jet, from the same image as shown in Fig[TJ displayed with a logarithmic 
response to enhance regions of low surface brightness. These VLA observations were carried out in 1985, when knot 
HST-1 was in a quiescent, pre-flare state, enabling the surrounding jet structure to be imaged more clearly. In 
addition to the bright knots HST-1, D, E, F, I and A, this image shows the twisted, helical filaments between the 
knots which LHE studied. 



knots, LHE found no significant differences between the filaments in the VLA and HST image. 
LHE interpreted these results as evidence for a twisted-helix structure, wrapped around the edge 
of the jet and emitting in radio and optical bands, between HST-1 and knot A. The observed radial 
difference in filament and knot location within the jet suggests that we can model the filaments 
without considering the mechanism responsible for the knots. 

Looking ahead to our KH analysis, the most important result from LHE is that the observed 
filament wavelength (projected on the sky) is not constant; it increases from Ap b ~ 2" ± 0.3" at 
HST-1 to A° b ~ 3" ± 0.3" at knot A. We also point out that the bright knots in the inner jet do not, 
in general, coincide with filament crossings. Although knot E appears to coincide with a filament 
crossing, Figure [2] shows that knots D, F, I and A are more complex; this is true in the optical as 
well as the radio images. 

2.2. Twisted Filaments and the KH Instability 

In what follows, we determine the implications for the jet and its immediate environment if 
the filaments are generated by KH instability. The KH instability can occur when velocity shear 
exists across a fluid boundary, such as the jet surface. It manifests itself in the form of pressure 
waves generated at or near the jet surface. The pressure waves can resemble the twisted filaments 
in the M87 jet. Small perturbations to the jet flow can be described mathematically in terms of a 
sum of Fourier components, referred to as "normal modes" (as in Section [4] and Appendix A). The 
two lowest-order twisting modes, helical and elliptical, describe associated distortions to the jet 
surface. Although higher -order modes i nitiall y grow more rapidly, the lower-order modes dominate 
at nonlinear levels, e.g., Hardee et al. (1997). Each mode involves both surface waves and body 



waves, which have different radial structures and move at different speeds. We will show (in Section 
14. 4p that surface waves dominate in the M87 situation. 

In this paper we argue that the elliptical mode is dominant in the inner jet, and gives rise 
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to the filaments described by LHE. We illustrate the pressure structure of the surface and first 
body elliptical modes in Figure O which shows a typical cross section and an integrated line of 
sight pressure image for each mode (jHarded 20001 ) . The elliptical surface mode generates two high- 
pressure filaments near the jet surface. The first body mode generates two high-pressure filaments 
located at about the radius midpoint, as well as a spiral, high-pressure pattern between these 
filaments and the jet surface. Because pressure can be used as a simple proxy for synchrotron 




Fig. 3. — Generic pressure cross sections (left) and line-of-sight pressure integrations (right) for elliptical surface 
(top) and elliptical first body (bottom) modes of a jet lying in the plane of the sky. Red indicates high and blue 
indicates low pressure. No relativistic or light-travel effects are included in the line-of-sight images. 

emissivity (see Section [6]), Figure [3] illustrates how the the surface and body elliptical modes might 
appear in a radio jet lying in the plane of the sky. Comparison of Figures [2] and [3l along with the 
double Gaussian fits from LHE, indicate that the filaments in the M 87 jet lie near to the jet surface 
and resemble the elliptical surface mode. 

We thus suggest that t he elliptical surface mode dominates in the inner M87 jet. Relativistic 
hydrodynamic simulations ( Hardee et al.l fl997l : lHardee Sz Hughes! 120031 ) suggest that the elliptical 
surface mode saturates without disrupting the jet, perhaps because this mode does not displace 
the jet's centerline. We note that the elliptical body mode grows more slowly than the surface 
mode in a linear analysis (see Section 14.41). and simulations show the body modes sa turate at lower 
amplitudes than do the surface modes ijHardee et al.ll200il ; lHardee &: Hughes! 120031 ). 

Our hypothesis that elliptical modes dominate in the inner M 87 jet disagrees with lBicknell &; Begelman 



( 19961 ) who suggested the bright knots in the M 87 jet are caused by the helical mode. In particular 
the inner M 87 jet bet ween HST-1 an d knot A (as in Figure [2]) does not resemble the morphology of 
the helical mode, e.g.. lHardee! ( 20001 ) . The helical mode displaces the jet centroid, causing a lateral 
oscillation in the jet; but the inner M 87 jet has very straight edges, with no sign of lateral displace- 
ment. In the absence of significant lateral oscillation, the helical surface mode could still appear as 
a single, high-pressure filament wrapped around the jet; this disagrees with the dual, double-helix 
filaments that exist in the real inner jet. By comparison, the elliptical mode does not displace the 
jet centroid, but it does create a pair of high-pressure twisted filaments w hich wrap around th e jet. 
Furthermore, no helical twisting has been detected at sub-parsec scales ( Acciari et al.l[2009l ): this 
suggests there is only a minimal helical perturbation to the jet inside HST-1. 

The apparent absence of a KH helical mode is perhaps surprising. Linear analysis (as in Section 
I4.4|) finds the helical mode will grow more slowly than higher-order modes in the M 87 jet, but sti ll 
rapi dly enough that one migh t expect it to be detectable. Simulations agree: lHardee et all (119971 ). 
also lHardee &: Hughes! ( 20031 ). find that the helical mode, if driven, grows slowly but surely, and 
does not saturate. Instead, it continues to grow until it reaches nonlinear levels and disrupts the 
jet. Thus, the lack of any detectable helical mode, either on sub-parsec scales or between HST-1 
and knot A, suggests that any initial helical-mode perturbation to the jet must be very small. 
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2.3. Jet Flow and Pattern Motions 

In order to use the observed filament wavelengths to determine conditions in the jet and its 
surroundings, we must determine the viewing angle, the jet flow speed and the speed of the wave 
pattern within the jet. To this end, we take advantage of the extensive data on proper motions of 
radio and optical features in the jet. Because we expect the KH instability to give jet flow through 
wave patterns, we identify faster proper motions with the jet flow, and slower proper motions with 
a pattern speed. The similarity of the optical and radio jet filaments, and the lack of evidence for 
a broad velocity shear layer, make it likely that we can use the fastest observed proper motions 
from both bands to reliably represent the jet speed. In what follows, we use proper motion data 
to develop continuous, steady-state models of the jet flow and pattern speed. By doing this, we 
implicitly assume that whatever causes the emission enhancement in the bright knots does not 
create any significant deviations from smooth flow within the jet. 

2.3.1. Viewing Angle 

To constrain the viewing angle, we look to the fastest observed proper motions. The fastest 
optical proper mot ion implies a superluminal speed /3 ob ~ 6.1 ± 0.6 at the posi tion of HST-1 
( Biretta et al. 19991 ). The most rapid radio proper motion implies /3 ob ~ 4.3 ±0.7 ( Cheung et al 



l2007|) at theposition of HST-1. Recalling that the maximum viewing angle is given by cos # n 
[(/3 ob ) 2 — l]/[(/3 ob ) 2 + 1], the highest optically determined superluminal speed requires a jet viewing 
angle 9 < 18.6°. On the other hand, the highest radio-determined superluminal speed requires only 
6 < 26.1°. We therefore work with two models, jets at 15° and 25° to the line of sight, taking these 
as representative of limits on the jet viewing angle and the jet speed at HST-1. 

2.3.2. Jet flow speed 

To constrain the flow speed of the jet plasma, we again look to the fastest observed proper 
motions. The fastest optical proper motions decline along the jet with superluminal speeds /3 ob ~ 
6.1±0.6,~ 5.11±0.66,~ 5.1±0.9, and ~ 3.9 ±0.8 found at the positions of HST-1, HST-2, knot D, 
and knot E, respectively ( Biretta et al. 1999). The fastest radio proper motions also decline along 



the jet with speeds of /3 ob ~ 4.3 ± 0.7 at HST-1 ([Cheung et al.1 120071 ) and (3 oh ~ 2.5 ± 0.3 at knot 



D dBiretta et alJll99fih. We us e the radio proper motion downstream of knot A, /3 ob ~ 1.25 ± 0.13 



at knot B ([Biretta et al.lll995l ). to set a lower limit on the flow speed at knot A. For HST-1 and 



knot D we observe flow through a more slowly moving knot. Thus, the fastest observed radio and 
optical proper motions associated with the knots suggest a flow speed which decreases along the 
jet between HST-1 and knot A. At a viewing angle of 15°, the four optical data points spanning 
the 6" from HST-1 to knot E require more than a 34% decrease in the Lorentz factor and provide 
the most robust evidence for a decreasing jet speed. 

We synthesize these data into two jet models. To represent the faster optical motions and a 
smaller viewing angle we work with a fast jet model with viewing angle 9 = 15° and intrinsic jet 
full opening angle 2ip = 0.0272 radian = 1.59°. The jet Lorentz factor is -y j0 = 7.50 (f3 oh ~ 6) 
at HST-1, 7j ~ 5.40 at knot D and jj ~ 4.15 at knot E. To represent the slower radio motions 
and a larger viewing angle we work with a slow jet model with 9 = 25°, full opening angle 2ip = 
0.0444 radian = 2.54°. In this model, the jet Lorentz factor is j j0 = 4.40 (/3 ob ~ 3.5) at HST-1, 
and 7j ~ 3.10 at knot D. Our fast and slow jet models are constrained by a lower limit to the jet 
Lorentz factor at knot A, jj > 1.90 for 9 = 15° and 7j > 1.70 for 9 = 25°. We return to these two 
jet models in Section [3l 

2.3.3. Observed pattern speed 

To constrain the pattern speed, we turn to slower, subluminal, optical or radio proper motions. 
Here our best evidence comes from the slower speeds associated with HST-1 or knot D, because these 
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two knots provide clear evidence for jet flow through more slowly moving structures (patterns), as 
might be expected for the KH instability. While knot D s hows a broad range o f motions, the average 
radio-determined speed for knot D is /3 ob ~ 0.41 ± 0.10 (iBiretta et al.lll995l ). The story at HST-1 



i s more complex. T he upstream edge of HS T-1 has been measur ed to move at /3 ob ~ 0.84 ±0.11 
Biretta et all (ll999T l. optical), at /3 ob < 0.25 (jCheung et al.1 (|2007l ). radio), and at /3 ob ~ 0.61 ±0.31 



Chang et al.1 ( 2010I ). radio). In addition to the known variability of HST-1, which makes it difficult 



to combine data taken at different epochs, radio data suggest the upstream edge of HST-1 may lie 
at the northern edge of the j et, possibly caused by an interaction between the jet and something 
in the external environment ( Cheung et al. 20071 ). We conclude that knot HST-1 is not the best 



choice for a pattern speed, and take the average radio-determined motion of knot D as the best 
available indicator of pattern speed. 

To determine the evolution of pattern speed along the jet, we assume sufficiently slow spatial 
growth so that the filament wavelength, X p , can be related to the pattern speed by X p ~ 2irv p /<jj. 
We assume that the wave frequency u is set at the filament source, somewhere upstream of HST-1, 
and remains constant along the jet. The observed 50% increase in \ p between HST-1 and knot A 
requires a 50% increase in v p over this distance. We tie this required increase to the average observed 
speed at knot D, /3 ob = 0.40. This choice requires "observed" pattern speeds of /3° b = 0.35 at HST-1 
and /3° h = 0.53 at knot A, which are consistent with the ranges in the data, namely /3 ob ~ 0.2-0.8, 
at HST-1 and /3 ob ~ 0.41 - 0.61 at knot A (IBiretta et al.lll995h . 



2.3.4- Intrinsic pattern properties 
Finally, we connect the observed pattern wavelength, A° b , to the intrinsic pattern wavelength, 
Ajf, by 

\ob _ ftp un _ \in /-t\ 

Xp ~ ^ " 1-/3^ cos A • (1) 

Figure U] shows the intrinsic pattern speeds and wavelengths from HST-1 to knot A accompanying 
the observed wavelengths and observed pattern speeds for the jet at 9 = 15° and at 9 = 25°. In 
general, the observed wavelength can be greater or less than the intrinsic wavelength depending 
on the intrinsic pattern speed, /3 p n , and viewing angle, 9. For our viewing angles the intrinsic 
wavelength and pattern speed are greater than the observed values, and show less change along 
the jet than the observed values. In particular, the intrinsic wavelength, A™, increases more slowly 

than the jet radius, Vj, so that AJf/Vj decreases along the jet. This has consequences for the KH 
instability analysis we present in Section 01 

3. Models of the Jet Flow Field 

In this section we develop our baseline jet models, which we need in order to analyze the behav- 
ior of the KH instability, and determine what conditions in the jet would accompany deceleration 
in accord with the observations. We consider a steady-state jet which conserves mass flux and is in 
pressure balance with an external cocoon. The cocoon acts as the ambient medium which controls 
the development of the KH instability. The full system conserves energy flux as the jet decelerates. 

We do not tackle the complex problem of what causes the jet to decelerate. Instead, we 
parameterize our jet models to describe the extent to which energy flux is conserved in the flow. 
We begin with models in which energy flux is conserved within the jet, so that the jet plasma heats 
up as it decelerates. We also explore jets in which a small fraction of the energy flux is transfered 
to the cocoon, as would be expected if the jet's interaction with the cocoon is responsible for the 
observed jet deceleration. Provided the energy lost from the jet is a small fraction of the total 
energy flux, both models give very similar results for the internal structure of the jet. 
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Fig. 4. — Observed and intrinsic wavelengths for the inner M87 jet. The observed wavelength Ap b (black dashed 
line) and error bars are from LHE; we fit a linear relation betw een these two points. The observed pattern speed f3° h 
(black dash-dot line), and error bar, from (|Biretta et al.|[l995l ). refer to knot D. Note that the pattern speed increase 
is set by the observed wavelength increase; thus we fit to only one data point. The intrinsic wavelength, Ap n (dashed 
lines), and pattern speeds, /3p n (dash-dot lines), are derived for our fast jet model (8 — 15°; green) and our slow jet 
model (8 = 25°; red), given the observed wavelength, A° b and the observed the pattern speed, /3p b . Knot positions 
and nomenclature are indicated by the letters above the lower axis. 



3.1. The Energy-Conserving Jet 

We first consider a jet which conserves mass flux and energy flux. Because the jet expands 
uniformly between HST-1 and knot A, we can assume flow along radial lines in spherical geometry. 
Conservation of mass flux can be written in terms of two locations along the jet, a fiducial point at 
axial position Zjo (with jet radius r^o) and another point at some general Zj (with jet radius rj): 

~.jf>.i V7 = "/(>/',(! ^o'jo » (2) 

where pj is the jet mass density; jj and f3j are the usual Lorentz factor and normalized flow speed. 
Conservation of energy flux can be written as 

IjWjfyr] = ^ W j0 p j0 fj , (3) 

where the total energy content of the plasma, the "relativistic" enthalpy (rest mass plus enthalpy), 
is described by 



J 

and Tj is the adiabatic index of the jet plasma. We also define the ratio of pressure to rest mass 

- w 



energy, 0, = Pj/pjc 2 , and the specific enthalpy, 



We note that all of the jet parameters - Jj(z), f3j(z), Pj(z), Pj(z), Tj(z), Xj( z ) an d @j(z) ~ are 
functions of position along the jet. We suppress this z-dependence to lighten the notation. 
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Values for Xj along an expanding jet containing relativistic plasma can vary from less than 
one to greater t han one. The ad i abatic index, als o changes with Xj hi tliis regime. We use the 
approximation ( Thompson Il986l : iDuncan et al.lli996l ) 



10.72 + 5pc 2 /P 
~ 8.04 + 3pc 2 /P ' 

and use this to connect Xj to 0, as 

Xi(©i) = ©i/(©i) = 6,(10.726, + 5) / (2.686,- + 2) . 



(6) 



(7) 



The last equality defines the order-unity function /(6), which has the expected limits, / (6) — > 2.5 
for 9 < 1, and /(6) ->■ 4 for 9 > 1. 

We now use Equation ([2|) in Equation (J3j) to rewrite energy flux conservation as 



7j 



1 + 



1 



Pj c 



7jo 



1 + 



jo 



jo 



1 



JO 



or 



(8) 



(9) 



H I 1 + Xj(@j)] = Ijo (1 + Xjo) , 
where Xjo = x(©jo)- This makes it apparent that a decelerating jet which conserves energy flux 
will heat up (7, < 7,0 =>• Xj > Xjo =^ ©j > @jo)) and vice versa. 

We want to solve Equation to describe the flow field in the jet. In order to do this we need 
to connect the enthalpy, Xjt to the flow field. Rather than modeling the full thermodynamics of 
the flow, we assume a polytropic dependence between the pressure and density, 



1 3 



Pjo) 



(10) 



for some (as yet unknown) index e ec . The subscript "ec" denotes "energy-conserving". For an 
expanding jet, e ec = Tj > 1 describes an accelerating adiabatic jet, e ec = 1 describes an isothermal 
jet at constant speed, and e cc < 1 describes a jet which heats up as it decelerates. 6ur polytropic 
assumption also gives 

£ec— 1 



_©_i_ 

0,o 



PL 

Pjo 



and 



Xj(©i) 
Xjo 



Ei 

PjO 



-1 



5(©, 



(11) 



(12) 



Here g(@j) = /(©j)//(@jo) is another order-unity function, which can be expressed in terms of 
(Pj/pjo) and xjo- 

3.2. Constraints on Pressure Decline Along the Jet 

Equations (fTT]) and (fT2j) show that the thermal state of the jet (described by x or 6) is uniquely 
related to 7, if the polytropic index e cc is known. In particular, e cc determines the pressure drop 
along the jet (from Equation [TUj). We therefore use the data to estimate s ec . Two sets of evidence 
are available: radio data which can be used to determine the minimum pressure in the jet and 
X-ray data which reveal the pressure in the local ISM. 
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Basic synchrotron theory can be used to find the minim u m jet pressure, P, 



Pacholczvk 



i m jet P 

(EOT; I 



mini 



consistent 
Burns et al.l (| 19781 ). OHC 



with the jet's observed luminosity and volume, e.g. 
calculated P m in from their radio observations, ignoring relativistic effects. They found that -P m in 
declines by a factor ~ 1.5 between HST-1 and knot D, and by a further factor ~ 2 between knots D 
and E. Because their data were obtained when HST-1 was in a quiescent state, their P m i n value for 
that knot should be typical of the underlying jet. The scatter in their P m i n estimates downstream 
of knot E suggests little additional decline to knot A, and thus an overall decline in Pm i n by a 
factor of ~ 3 between HST-1 and knot A. Similar calculations carried out bv lBiretta et al. ( 19911 ). 
including optical and X-ray data, obtained similar results. 

Stawarz et al.l ( 20061 ) revisited P m i n in the jet, including relativistic effects via a beaming 
correction P m ; n oc <5~ 10 / 7 (jStawarz et al.l [20031 ). where the Doppler factor is 5 = [7(1 — /Jcos^)]" 1 . 
In Section \2. 31 we argued that typical jet Lorentz factors ~ 4.4—7.5 at HST-1, slowing to ~ 3.1—5.4 
at knot D. These values, along with their respective viewing angles 25° and 15°, give Doppler 
factors 5 ~ 1.9 — 3.1 at HST-1 and 5 ~ 2.2 — 3.6 at knot D. (The slightly larger values at knot D 
are a consequence of the /3cos9 dependence in 5.) This result suggests an additional decline in 
-Pmin) relative to the uncorrected estimate, of ~ 20 — 25% from HST-1 to knot A. Overall, these 
calculations suggest that P m i n declines by a factor ~ 3 — 4 between HST-1 and knot A. There is, 
of course, no fundamental reason why the pressure of the jet plasma should be equal to P m i n - In 
particular, if the jet is only weakly magnetized, we expect Pj > P m \ n . Nonetheless, we might guess 
that a decline in P m i n approximately reflects the decline in Pj. 

The evolving jet and cocoon need not be in pressure balance with the ambient ISM and mini- 
mum pressure estimates (Section I5.2f) indicate an overpressured jet. Even so, it is still interesting 
to consider the pressure drop in the ISM along the jet. On ~ 10 — 100 kpc sca les, the ISM is fairly 
smoot hly distributed, and its pressure declines smoothly with radius, e.g., iNulsen fe; Bohringerl 
(j 19951 ). However, the ISM in the inner part of the galaxy is complicated (as in Section [2]), 
and the ISM pressure pr obably fluctuates substantially about any large-scale trend. Nonethe- 
less, [Stawaxz^et^aL ( 20061 ) estimated a factor ~ 3 decline in the ISM pressure between HST-1 and 
knot A which is interestingly similar to the factor ~ 3 — 4 drop in P m i n from the radio data. 

Based on this evidence, we proceed by modeling a jet whose pressure declines by a factor ~ 3— 4 
between HST-1 and knot A. The value of e e c connects the pressure decline to the jet flow field. 
Using Equation (110p . we find e cc = 0.325 ± 0.025 provides a pressure decline ~ 3.5 ± 0.5 between 
HST-1 and knot A. This result is nearly independent of the initial jet speed and viewing angle used 
to fit the observed jet deceleration. In what follows we set e cc = 0.325 for our energy-conserving 
jet models. Note, we do not specify an absolute value for the pressure at HST-1; our KH analysis 
only needs the pressure change along the jet. We discuss pressure scaling further in Section 



3.3. Structure of the Energy- Conserving Jet 

Having chosen the polytropic index, we can build jet models which conserve energy flux and 
mass flux. In Section [3.11 we showed that the two parameters, Xjo an d e e o totally specify the flow 
field in the jet, if the Lorentz factor, 7^0, is known at some initial point rjQ. To solve for the flow 
field, we combine equations and (fT^j) . as 



7jo(l + Xjo) =7j 



(13) 



This equation has only one unknown, jj, the jet Lorentz factor at radius rj = ipZj (recall that tp is 
held constant and g(Qj) is an order-unity function which can be written in terms of the jet speed, 
radius, Xjo an d £ e c)- I n our solution, we begin at HST-1, choosing e cc = 0.325, and a specific set of 
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initial conditions (rjo, Tjo, Xjo)- We solve Equation (fT3|) numerically, to determine Jj(z) at several 
positions along the jet, out to knot A (where rj ~ 12rjo)- We compare the results to the values of 
/3? b at HST-1, knot D, knot E and the lower limit to /3? b at knot A, then adjust y-io to iterate on 
Equation (fT3|) as needed, to give the best overall fit to the data for our chosen e cc . Finally, we use 
Equations (jlOp or (jlip to find the evolution of the internal energy along the jet. 

In Figure [5] we show our best fits to fast and slow energy-conserving jets. Our fast jet with 
9 = 15° starts at HST-1 with 7^0 = 7.50 and Xjo = 0.130. By knot A it has decelerated to jj = 2.66 
and heated to Xj = 2.189. This model easily fits the optical superluminal speeds and significantly 
exceeds the lower limit at knot A. Our slow jet with 6 = 25° starts at HST-1 with 7^0 = 4.40 and 
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Fig. 5. — Structure of the energy- conserving jet, for our fast jet (green lines) and slow jet (red lines) models. Left: 
the observed jet speed speed, /3° b , (dashed lines) and intrinsic jet Lorentz factor, jj, (solid lines). The optical speeds, 
/? ob , at HST-1, knot D and knot E are indicated by filled squares and radio speeds, f3 oh , at HST-1 and knot D are 
indicated by filled circles, both with error bars. The filled circle and arrow at knot A indicate the lower limit to 
/? ob set by the fastest observed radio superluminal speed at knot B. Right: the specific enthalphy, Xj> normalized 
temperature Qj = Pj/pjC 2 , pressure Pj/Pjo and density Pj/pjo, from solutions to Equation (|13[1 with e oc = 0.325. 



Xjo = 0.095. By knot A it has decelerated to 7j = 1.82 and heated to Xj = 1-643. This model is just 
barely acceptable; it lies just outside the error bars associated with the fastest radio superluminal 
speeds at HST-1 and knot D, and just above the lower limit at knot A. We therefore treat our 
slow-jet model as a lower limit to the jet speed and upper limit to the viewing angle. In both fast 
and slow models the enthalpy rises from ~ 0.1 at HST-1 to ~ 2.0 at knot A; the jet plasma heats as 
it absorbs the lost kinetic energy. The jet plasma is hot but subrelativistic, P ~ 0.04/3C 2 , at HST-1 
and heats by a factor ~ 15 to P ~ 0.6pc 2 at knot A. We do not derive numerical values for the 
pressure or density, but only the ratios of each quantity to its value at HST-1; these are also shown 
in Figure [5j Although the jet decelerates, its density still drops due to the expansion, by a factor 
~ 50 between HST-1 and knot A. The jet pressure is, of course, constrained by our assumptions to 
drop by a factor ~ 3.5 over the same range. 

Solutions to Equation (|T3|) are sensitive to the starting value of the enthalpy, Xjo- Flows in 
which the jet plasma is internally relativistic at HST-1 (xjo S> 1) decelerate very rapidly and violate 
the lower limit on jet speed at knot A. Flows which are very cool at HST-1 (xjo ^ 1) decelerate 
too slowly to match the proper motions through knots D and E. Specifically, we find an upper 
limit to the starting enthalpy, Xjo ax = 0.74 by considering a maximum allowed deceleration from 

7j0 = 10.5 at HST-1 to jj = 1.8 at knot A (using the upper error bar limit, /3 ob < 6.7 at HST-1 
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and the lower error bar limit /3 ob > 1.12 at knot A). We also find a lower limit of x^ in = 0.032, by 
considering a minimum allowed deceleration from 7^0 = 6.3 at HST-1 to 7j = 4.7 at knot E (using 
the lower error bar limit, (3 oh > 5.5, at HST-1 and the upper error bar limit, /3 ob < 4.5, at knot E). 
The minimum deceleration case gives jj = 3.3 at knot A. 

3.4. The Energy-Losing Jet 

In our energy-conserving jet, all of the excess kinetic energy lost as the jet decelerates is con- 
verted into internal energy within the jet. Because the jet probably decelerates through interactions 
with the surrounding medium, some of that excess kinetic energy may well go into the surroundings. 
This could happen, for instance, as a consequence of shocks driven into the cocoon by the unstable 
KH modes. We want to explore how such energy loss affects the KH instability. However, building 
a complete model of jet and ISM energetics is beyond the scope of this paper. Instead, we invent 
a toy model to describe the effect of the energy transferred from the jet to the cocoon. 

We know that the pressure in an energy- losing jet must drop more rapidly, relative to the local 
density, than in an energy-conserving jet. However, the "observed" jet pressure must still drop by 
a factor ~ 3.5, to match our arguments in Section 13.21 We parameterize this as follows. We first 
assume that the underlying jet flow satisfies equations ([2]) and (|3|), and the pressure also obeys a 
polytropic relation, but now with an "energy-losing" polytropic index, e c i < e ec = 0.325. We solve 
Equation (|13p with this new e e \, and match the solution to the observed velocities. We adopt the 
same Lorentz factors at HST-1 (7.50, fast jet, and 4.40, slow jet), and at knot A (2.66, fast jet and 
1.82, slow jet), as for an energy-conserving jet, and we again adjust Xjo to match these data. The 
values of e e i and Xjo consistent with the observed deceleration are interrelated; smaller values of 
Xjo accompany smaller values of e e i. At this point, to ensure that the "true" jet pressure drops 
more rapidly than e e i describes, we assume the "true" jet pressure is related to the jet density 
via Equation (110p . but now with the energy-conserving e ec = 0.325. Thus, we parameterize the 
mismatch between the "observed" pressure decline, and an energy-conserving pressure decline, by 
e ec — e e \. We return to these models in Section [4.21 where we will show that energy-losing models 
consistent with the data must have e e i > 0.250, and that models with 0.250 < e e i < 0.325 differ 
only slightly from the energy-conserving models. 

4. KH Instabilities in the Jet-Cocoon System 

In this section we determine under what conditions the twisted filaments seen in the M87 
jet can be caused by the KH instability. We begin with a brief overview of the KH instability, 
including some useful analytic approximations; we store more details in the Appendix. We analyse 
our energy-conserving and energy-losing jet models, finding what conditions must exist in the 
cocoon if the KH-generated filaments in our model jets are to match the ones in the real jet. We 
show that the KH modes propagate supersonically in the jet, and verify that they grow rapidly 
enough to saturate close to the origin of the jet. 

KH stability analysis describes the growth of perturbations to the jet using solutions of the 
dispersion relation, Equation (|A2p in the Appendix. This equation can be solved to determine the 
wavelength, A^, and wave speed, v w , associated with a frequency, u w . The dispersion relation 
solution describes a 180° rotation of the elliptical cross section distortion, corresponding to a 
180° rotation of a single filament. Thus, the elliptical wave frequency and wavelength described 
in solutions to the dispersion relation are related to the observed frequency and wavelength by 
oj w = 2uj p and A^ = A p /2. 

Numerically, we rely on the intrinsic pattern wavelength, A p , and speed, v p , derived from the 
data (as in Section [2]), to determine the associated pattern frequency, uj p = 2ttv p /X p . Specifically, 
from Figure [H the intrinsic pattern speeds at HST-1 are v p ~ 0.60c for the fast jet, and v p ~ 0.48c 
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for the slow jet; the intrinsic wavelengths are A p ~ 3.3" for the fast jet and ~ 2.7" for the slow 
jet. From these we estimate u p ~ 1.4 x 10~ 10 rad s _1 , corresponding to a periodicity of ~ 1440 yr. 
We assume uj p is set at the filament source, somewhere upstream of HST-1, and remains constant 
along the jet. 

We assume the jet is in pressure balance with the cocoon, as required by our KH analysis. 
Solutions to the dispersion relation (see details in the Appendix) depend on the flow speed and 
thermal state of the jet and the cocoon. In practice, at a given location along the jet, we use the jet 
parameters from our models in Section [3] and iterate on the cocoon parameters until we find a good 
solution. Specifically, we begin with an initial guess for the cocoon sound speed, solve Equation 
(|A2j) numerically for the complex wavenumber and compute the associated wavelength A(w) and 
wave speed, v{ui) as a function of oj over a broad frequency range. We iterate on a c until we obtain 
a value of v w = v p at the wavelength X w = X p /2 that is within 2% of the observed values. This 
typically requires five to ten iterations. Note that we are matching computed wavelengths and wave 
speeds to observed intrinsic wavelengths and wave speeds, and the frequency will be ui w = ui p /2. 
We do this at seven locations along the jet: HST-1, 2", knot D, knot E, knot F, knot I and knot A. 

4.1. Analytic Approximations Close to Resonance 

KH instabilites are characterized by a resonant frequency and wavelength at which the growth 
rate is a maximum. In general we expect perturbations to the jet to produce structures close to 
this resonance. In Section 14.21 we find that the elliptical surface mode in our model jets is close to 
resonance. Had we found the alternative - that the observed filaments were not close to resonance 
- our identification of the filaments with KH instability would be called into question. 

While the KH dispersion relation must be solved numerically, as we do in Sections 14.21 and 
14.31 analytic expressions for the resonant conditions provide insight into the numerical results. The 
resonant solutions, given in Equations (|A4|) and (|A3|) in the Appendix, depend on conditions in 
the jet and the cocoon. These expressions can be simplified for conditions found in our numerical 
solutions, namely, u c <C Uj ~ c. In this limit the resonant wave (pattern) speed becomes 

< ~ (14) 



Uj JjClc + CLj 

Here, Uj = f3jc is the jet speed, 7? = (1 — , and dj and a c are the sound speeds in the jet and 
cocoon. The sound speed depends on the thermal state of the plasma: 

^ = ^ = (r " 1} JTT7) = i + re°r-i) ' (15) 

In Equation (j!4j) we have also used the fact that j a ~ 1 (because a < 0.5c). From Equation (|14p 
we verify that resonant waves move more slowly than the underlying flow: ?;* < Uj. We also see 
that a change in the intrinsic wave speed, v^, requires a change in the cocoon sound speed, a c , 
which depends on the jet model through jj and cij. The resonant wave frequency becomes 



-1/2 



(16) 



a c 4 

where r,- is the local jet radius. We can combine Equations (|14p and (|16p with = 2-rrv^/u}* to get 

A* 8 Uj [1 - (a c /v*J 2 ^ 



5a c [1 + (aj/jja c )] 



(17) 
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Equation ()17|) shows that the cocoon thermal state - described here by the sound speed - is the 
only important parameter (given u c <C a c ). If A* jvj decreases, Uj/a c must decrease. But, because 
uj ~ c, a c must rise in order to account for the filaments in terms of a KH instability. 

4.2. Limiting-Case Jet and Cocoon Models 

In addition to energy-conserving jets (Section 13. 3|) . we want to analyze the energy-losing jets 
(Section I3.4j) . In the latter case the cocoon as well as the jet is involved in conservation of energy 
flux, so we replace Equation ([3|) by 

j 2 c W c f3 c (r 2 c - rj) + -'jWyh'j = + TcO^cO&o^o - r%) (18) 

where the subscript "c" labels properties in the cocoon. For the energy-losing jet models, we self- 
consistently compute f3 c from Equation (|18p . as part of the iteration to solve for a c . We assume 
that the observed jet carries essentially all the energy flux at HST-1: 7co^ / co/?co( r cO ~~ r jo) ^ 
Ijo^joPjOfjo- We find that the cocoon flow speed is dynamically insignificant, /3 C <C f3 ac , provided 
the cocoon radius r c > 4rj; in specific models we assume r c = Wrj. 

We can only find physical solutions to cocoon conditions in the energy-losing jets if 0.250 < 
e e \ < 0.325. At smaller values of e e \ the cocoon sound speed necessary to match the filament 
properties goes to zero somewhere between HST-1 and knot E, requiring the cocoon density to 
become infinite in order to maintain pressure balance. The lower limit on e e i corresponds to an 
upper limit on energy lost from the jet to the cocoon, ~ 19% for fast jets and ~ 16% for slow jets. 
In this section we take e i = 0.250 for our energy-losing jets, to determine the effects of the largest 
possible energy loss. 

We show the results of our KH analysis in Table 1 and in Figure [6j Table 1 gives values at 
the fiducial point, HST-1, for our four limiting-case jet models (fast and slow, energy-conserving 
and energy- losing). The thermal state of the jet is described via the enthalpy, Xjo, the adiabatic 
index, Tjo, and the sound speed f3 a j,o- The thermal state of the cocoon is given via the same three 
quantities. Table 1 also contains the wave frequency, u) w , and wavelength, X w , that fits the observed 

Table 1: Jet, Cocoon & Elliptical Mode Parameters at HST-1 



Model 


z(pc) 


Xjo 


r i0 


Ay',0 


XcO 


r c o 


Pac,0 


uj w rj/uj 


Xyj/Vj 




- 1 - ec 
Fel 


312 
312 


0.130 
0.092 


1.627 
1.637 


0.268 
0.232 


0.0064 
0.0038 


1.664 
1.665 


0.065 
0.050 


0.119 
0.119 


31.5 
31.5 


31.5 
38.6 


Sec 
Scl 


192 
192 


0.095 
0.068 


1.637 
1.644 


0.235 
0.203 


0.0079 
0.0046 


1.664 
1.665 


0.072 
0.055 


0.122 
0.122 


25.5 
25.5 


23.4 
28.5 



intrinsic wavelength and wave speed along with the resonant wavelength, A^, of the elliptical KH 
mode computed from the dynamical and thermal state of the jet and cocoon. Quantities are scaled 
to the jet radius, rj, and speed, Uj, at HST-1. We see that the intrinsic filament wavelengths at 
HST-1, if produced by the elliptical surface mode, are within 20% of the resonant wavelength for 
that mode. This is in accord with our expectation that perturbations will produce structures close 
to the resonant wavelength. 

In the upper panels of Figure [6] we show the density, pressure and internal energy for both fast 
and slow versions of both jet models, interpolated with a seventh-order polynomial between the 
numerical solutions at the seven locations along the jet. The middle panels of Figure [6] show the 
cocoon density and temperature structure required to explain the filaments in the jet in terms of 
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the KH instability. Because we assume pressure balance, the pressure decline in the upper panels 
of Figure [6] also applies to the cocoon. The lower panels of Figure [6] highlight the internal energy in 
the jet and cocoon, via the sound speed in each plasma. In both fast and slow models, the cocoon 
starts cooler than the jet, but heats up more rapidly, so that a c ~ aj by knot A. 
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Fig. 6. — Jet and cocoon parameters for the fast jet models (left) and slow jet models (right), as a function of 
intrinsic distance along the jet. Energy-conserving jets are shown in blue; maximum energy-losing jets in magenta. 
The upper panels show the jet temperature, relative pressure and relative number density. The middle panels show 
the cocoon temperature and relative number density. The bottom panels shows sound speeds in the jet, aj/c, and 
cocoon, a c /c. 



One key result is the thermal state of the cocoon required in order for the observed filaments 
to be caused by the KH instability. At HST-1, the cocoon must be cooler than the jet, but still 
quite warm; XcO ~ 0.004 — 0.008, or P c q/ Pcoc 2 ~ 0.002 — 0.003. Because the jet and cocoon are in 
pressure balance, the cocoon must be denser than the jet at HST-1: p c o/PjO — Xjo/XcO ~ 12 — 24. 
By knot A the internal energy of the cocoon must be nearly relativistic, and comparable that of 
the jet: \cA ~ 1-3 — 2.4, or P c a/PcAC 2 ~ 0.4 — 0.7. Pressure balance requires p c A ~ PjA at knot A. 

We emphasize that the increase in cocoon temperature along the jet is not a result of energy 
deposition from the decelerating jet. Rather, the temperature increase is required to explain the 
observed filament wavelength, which increases less rapidly than the jet radius (e.g., Figure 0]). As 
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discussed in Section I4TT] the decrease of X w /rj along the flow requires Uj/a c to decrease as well, 
and thus a c must increase because Uj — c throughout the jet. We do not try to model the thermal 
state of the cocoon; such modeling is beyond the scope of this paper. We simply note that such a 
hot cocoon must exist if the observed filaments are caused by the KH instability. 

4.3. Intermediate-Case Jet and Cocoon Models 

Figure [6] shows only small differences between the two limiting cases of energy-conserving 
and energy-losing jets. To continue, therefore, we develop an intermediate case, £i nt = 0.300, 
corresponding to an energy loss of 7% for the fast jet and 6% for the slow jet. We develop fast and 
slow flow models for this intermediate-case jet, as in Section [3. 4| and use these models for specific 
calculations in the rest of the paper. We carry out the filament-matching KH analysis, as in Section 
14. 2\ we give the details of these models in Table 2. We again find the filaments are close to resonance 

Table 2: Jet, Cocoon and Elliptical Mode Parameters for the Intermediate-Case Jet 
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z(pc) 
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iHST 
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0.991 
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1.000 
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0.006 


1.000 


0.061 


0.119 


31.5 


33.5 


F 


m" 
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0.988 


0.285 


0.289 


0.359 


0.010 


0.386 


0.080 


0.239 


16.0 


22.2 


F 


ikntD 


1090 


0.982 


0.567 


0.120 


0.439 


0.019 


0.150 


0.110 


0.420 


9.41 


14.8 


F 


ikntE 


1880 


0.970 


1.048 


0.052 


0.501 


0.078 


0.028 


0.215 


0.730 


5.70 


8.19 


F 


ikntF 


2660 


0.954 


1.517 


0.033 


0.530 


0.279 


0.007 


0.360 


1.05 


4.16 


4.38 


F ( 


ikntl 


3280 


0.939 


1.878 


0.025 


0.543 


0.642 


0.003 


0.460 


1.32 


3.44 


2.48 


F ( 


ikntA 


3750 


0.927 


2.140 


0.021 


0.550 


1.291 


0.002 


0.525 


1.53 


3.06 




S t 


mST 


192 


0.974 


0.085 


1.000 


0.224 


0.007 


1.000 


0.068 


0.122 


25.5 


24.9 


s - 


n" 


383 


0.967 


0.206 


0.282 


0.323 


0.012 


0.371 


0.090 


0.246 


13.1 


16.0 


St 


IkntD 


670 


0.954 


0.405 


0.110 


0.406 


0.024 


0.138 


0.125 


0.436 


7.73 


10.5 


S®kntE 


1150 


0.926 


0.737 


0.049 


0.475 


0.120 


0.022 


0.260 


0.770 


4.75 


5.23 


s, 


IkntF 


1630 


0.892 


1.052 


0.030 


0.509 


0.355 


0.007 


0.390 


1.13 


3.50 


1.83 


s, 


ikntl 


2010 


0.861 


1.288 


0.023 


0.525 


0.858 


0.003 


0.490 


1.45 


2.93 




s, 


'tkntA 


2300 


0.836 


1.454 


0.020 


0.534 


2.381 


0.001 


0.560 


1.71 


2.62 





over most of the jet. However, Table 2 does not give a reso nant wave length at knot A for either 
case, or at knot I for the slow jet case. Stability analysis ( Gill 19651 ) shows that the resonance 



which leads to a maximum in the growth rate does not exist below some supersonic limit which 
depends on the jet Lorentz factor and sound speeds. For example, Equations (|A3p or ([To]) show a 
rapid rise in the resonant frequency as o c — > < Uj (noting that 7j(7a )C a c ) ~ (TajOj)). Because 
uj* — > oo in this situation, the growth rate continues to increase in accord with the low-frequency 
approximation, Equation (jA5|) . at all frequencies. 

We can use this intermediate-case jet to show that the elliptical mode speed is supersonic in 
both cocoon and jet at HST-1 and is weakly supersonic in both the cocoon and the jet at knot I. 
That is, the jet plasma moves supersonically through the high pressure filaments associated with 
the elliptical mode, and the filaments themselves move supersonically relative to the cocoon plasma. 
We demonstrate this in Table 3 for fast and slow versions of the intermediate jet where we give the 
relativistic Mach numbers for the wave speed relative to the cocoon, M w<c = ('Jw,cPw,c)/{ r Ya,cl3a,c)} 
and for jet flow through the wave in the jet reference frame, Mj W = {"fj.wPj.w) / {la,jfia,j)- Thus, 



Table 3: Wave Speeds and Mach Numbers for the Intermediate-Case Jet 



Model 


Pj 


Pw,c 


Pj,w 


Paj 


Pa,c 


Af i>w 


M w , c 


F@HST 


0.991 


0.593 


0.965 


0.255 


0.061 


14.0 


12.1 


F@D 


0.982 


0.620 


0.925 


0.439 


0.110 


5.0 


7.2 


F@E 


0.970 


0.644 


0.869 


0.501 


0.215 


3.0 


3.8 


F@F 


0.954 


0.665 


0.791 


0.530 


0.360 


2.1 


2.3 


F@j 


0.939 


0.681 


0.716 


0.543 


0.543 


1.6 


1.8 


S@HST 


0.974 


0.480 


0.928 


0.224 


0.068 


10.8 


9.0 


S@D 


0.954 


0.509 


0.865 


0.406 


0.125 


3.9 


4.7 


S@£ 


0.926 


0.536 


0.774 


0.475 


0.260 


2.3 


2.3 


S@F 


0.892 


0.561 


0.663 


0.509 


0.390 


1.5 


1.6 


S@I 


0.861 


0.579 


0.562 


0.525 


0.490 


1.1 


1.3 



the KH wave mode can produce supersonic shocks in both the cocoon and the jet at HST-1, and 
weakly supersonic shocks at knot I or knot F. This supports our conjecture that the KH instability 
in the jet decelerates and heats the jet, and dissipates some energy in the cocoon, by means of 
shocks generated by the filaments as they reach large amplitudes. 

4.4. Normal Mode Growth Rates 

In the previous sections we determined conditions that must exist in the cocoon if the elliptical 
surface mode of the KH instability is responsible for the filaments in the M87 jet. In this section 
we verify that the growth rate of this mode is sufficiently rapid to saturate along the jet. We 
also investigate whether other modes, in harmonic resonance with the elliptical mode, might have 
comparable growth rates. We work with the two intermediate-case jet models from Section \4. 31 We 
use the full dispersion relation to compute the spatial growth length, £(z) = |/cj(,z)| _1 , where ki, the 
spatial growth rate, is the imaginary part of the wavenumber associated with the elliptical mode 
frequency, ui w = 2ttv w /X w . We use continuously variable conditions in the jet and cocoon found by 
interpolation between the seven locations shown in Table 2. We also calculate the spatial growth 
rates of the helical, triangular and rectangular modes, at the harmonically resonant frequencies 
uh = c*>u>/2, oje = w w , ojt = 3w^/2 and ujr = 2oj w . These harmonic frequencies produce filaments 
with wavelength and wave speed comparable, but not identical, to the elliptical mode. We chose 
these frequencies based on numerical simulations w hich show that othe r normal modes are excited 
in frequency resonance with the dominant mode ( Hardee et al.l [T997I ). which we assume is the 



elliptical surface mode. Because we have found that the elliptical surface mode is operating near to 
resonance along the jet, the assumption that other surface modes operate at the specified frequency 
resonances guarantees that these surface modes are close to their fastest growing frequencies along 
the jet, i.e., w* ~ nto* H where n = 1,2, ... is the helical, elliptical, etc. mode number. 

In Figure [7] we plot the spatial growth rate, \ki\ normalized by rj = tpZj, as a function of 
position along the jet. We find that each surface mode grows more rapidly than the accompanying 
body mode along the jet. The spatial growth rate of the surface modes when scaled to r,o at HST-1 
decreases slightly along the jet. The triangular mode growth rate is not shown in this figure, for 
reasons of clarity; it lies between the rectangular and elliptical mode growth rate. All four modes, 
however, are used in the pseudo-synchrotron image calculated in Section [6l Figure [7] also shows 
that surface modes will develop rapidly along the M87 jet; growth lengths for the elliptical surface 
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Fig. 7. — Spatial growth rates, \ki\rj, of harmonically resonant rectangular (R), elliptical (E) and helical (H) modes, 
as a function of distance along the jet. Surface (s) modes are shown by solid lines; body (b) modes are shown by 
dashed lines. The left panel shows the fast jet and the right panel shows the slow jet, both for the intermediate-jet 
model. 



mode range from £ > lOr^o at HST-1 to £ > rj = VItaq at knot A. We therefore speculate that 
surface modes can grow to significant amplitudes. 

Although we find that higher order surface modes initially grow more rapidly than the elliptical 
mode, we argue that the elliptical mode dominates in the actual M87 jet. The real jet must be 
bounded by an outer velocity shear layer. Linear analysis shows that hi gher-order modes c an be 
significantly affected, or even stabilized, by the presence of a shear layer ( Ferrari et al. 19821 ). but 



the helical and el liptical surface modes remain unstable even in the presence of significant shear 
( Birkinshawl Il99ll ) . In addition, simulations show that higher-order modes initially grow rapidly, 



but then saturate and are over whelmed by more slowly growing, lower-order normal modes that 
saturate at a larger amplitude ( Hardee et al.|[l997l ). The spatial growth rate of the helical surface 



mode is typically half that of the elliptical surface mode. This is enough for the helical mode to 
overwhelm a saturated elliptical mode, if it is excited at all. Thus, any initial helical perturbation 
to the M 87 jet must have been very tiny. 

We find that frequency-resonant body modes are also unstable in the inner region of the jet 
(Figure [7]) but they grow more slowly than surface modes. The very low growth rates of the body 
modes around HST-1, and also in the outer part of the jet, occur because the body modes are near 
their minimum unstable frequency in these regions. Complete stabilization of the body modes at 
all frequencies occurs between knots I and A for the fast jet model and between kn ots F and I fo r 
the slow jet model, because the jet is not sufficiently supersonic in those regions (|Hardedl2007l ). 



Based on this we argue that body modes do not develop to significant amplitudes on the M87 jet. 

5. An Example: Possible Jet and Cocoon Compositions 

In this section we explore an example of the possible internal state of the jet and the cocoon. 
Recall that our models do not derive the density or pressure of either system, but only the specific 
internal energy of each. Without definite information on the composition of the jet or cocoon 
plasmas, we cannot go any further. However, we can gain insight by considering an example in 
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which both the jet and cocoon contain a mix of relativistic and thermal plasmas, both scaled to 
the likely pressure in the M87 jet. 

5.1. Composition: Mixed Plasmas 

Let us consider a mixed plasma, which could describe the jet or the cocoon. Let it contain 
"thermal" ions, with some temperature T so that Pth = n^kT. Let it also contain relativistic 
leptons, with P rc \ = j e gn Te \m e c 2 , where j e s S> 1 describes the effective temperature of the leptons, 
normalized to m e c 2 . Note this differs from the bulk flow Lorentz factor of the plasma, which is 
much smaller. The leptons inside the jet are observed by means of their synchrotron radiation. The 
ions could provide charge balance with the relativistic leptons, or could represent a separate, ion- 
electron population. For the cocoon, the thermal component could have its origins in the galactic 
ISM, and the relativistic component could have its origins in the jet plasma, now mixed with the 
inner ISM and visible as the inner radio halo. 

Our models show the plasma obeys P = Qpc 2 , with O < 1. From this, 

P = P rci + P th = n vc a cS m e c 2 + n th kT = n rc] Qm e c 2 + n th @m p c 2 . (19) 

We guess that the relativistic leptons are well separated in energy from the thermal plasma. For 
instance, we know the jet spectrum continues to rise down to ~ 300 MHz (lOwen et al.ll2000h . If 
this were a low-frequency cutoff to the jet spectrum, and if the jet magnetic field £ 100/iG (as in 
Section [5T2]) the lepton pressure would be dominated by particles at 7 e g- > 300, so that 7 e g- S> 0. 

We normalize the ion temperature as £ = kT/m p c 2 , with no constraints on T except that 
( < is required to solve Equation (fT9j) . Solving that equation gives 



— ~ — — ; ~ p - (20) 

Ah C «th 7cff rn e 

At HST-1, where 0j ~ 0.04, and C ~ 0.002, the thermal plasma dominates the number density 
in both jet and cocoon (for any 7 e g >, 100). The pressure ratio depends on the unconstrained 
factor C; but unless the thermal plasma is much hotter than the T ~ 10 K typical of the galactic 
ISM, both jet and cocoon pressures are dominated by the relativistic components. At knot A, 
where Qj ~ C ~ 0.5, the detailed numbers may change but the situation is probably the same: 
the densities of the cocoon and jet plasmas are dominated by the thermal component, but their 
pressures are dominated by the relativistic component. 

5.2. Pressure, Density and Magnetic Field 

Up to this point our modeling has not depended on absolute values for the density or the 
pressure but only on the ratio of these two quantities. We can estimate the jet or cocoon densities 
by choosing one more parameter, namely, the jet pressure. From Equation (|19p we find that 

0-C P P 7cff ~ f0U 
nrd = — 7 tt 2 > nth = 2 7 ( 21 > 

(7cff - Qrn e c z Om p c z 7 cff - ( 

To estimate the pressure at our fiducial point, HST-1, we start with the minimum-pressure estimate 
from OHC, P min ~ 2 x 10" 8 dyn cm 2 . This value will be reduced a bit by the 5 10 / 7 Doppler 
correction, but the true jet pressure must exceed P m i n . We thus scale the pressure at HST-1 as 
Pjo = 10~ 8 i-b,8 dyn cm~ 2 . We guess that the relativistic leptons at HST-1 have 7j )C ff ~ 300, as 
above. We simply scale the relativistic component of the cocoon to 7 c ff iC = 1007 Cj ioo, for lack of 
any better information about relativistic leptons in the cocoon. 
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Putting in numbers, we start at HST-1, with Qj ~ 0.04 for the jet; this gives rij jTe \ — 4 x 
10~ 5 P ,8 cm" 3 , and n jjth ~ 2 x 10" 4 P 0j8 cm" 3 . For the cocoon at HST-1, with 9 C ~ .002, we 
get n ct h ~ 3 x 10~ 3 Po,8 cm" 3 and n Cjre i ~ 1 x 10~ 4 -Po,8/7c.ioo cm" 3 . At knot A, our models say 
the pressure drops by a factor 3.5 relative to that at HST-1. Both the jet and the cocoon have 
C ~ 0.5. If 7j jC ff °c @j holds for the jet, as above, we have Jj >e ft ~ 3800 at knot A. With 



these numbers, the jet densities are nj^ re \ ~ 1 x 10' 

-6n „~,-3 



Po,8 cm 3 > an d rij 



l j,th ~ 4 x 10~ b P ,8 cm" 3 . For 

the cocoon, we have n C) th ~ 4 x 10 _B Po,8 cm -3 , and n c ^ c \ ~ 4 x 10~ 5 Po,8/7c,ioo cm -3 at knot A. 
Thus, our assumed cocoon is not only hot, but also tenuous and over pressured, when compared 
to the n t h ~ 0.1 cm -3 and T ~ 10 7 K typical of the X-ray loud ISM in the inner few kpc of the 
galaxy. We discuss this cocoon further in Section fe.ll 

Our models assume the plasma is only weakly magnetized, Pb = B 2 /8n <C Pj- Numerically, 

1/2 

this gives B <C 500P 8 ' >G. Because Pj ~ P Tci in our example, our models must also have Pg <C P rc \- 
At first glance, this seems to contradict the results of Stawarz et al. (2005), who used upper limits 
on inverse Compton emission, together with minimum-energy analysis of the radio power, to argue 
PB^Prei m knot A. We note, however, that different assumptions in such modeling (compare, e.g., 
OHC), as well as smaller assumed viewing angles (as in Section T2.3p can accomodate a weaker field 
throughout the jet. In addition, Stawarz et al. (2005) assumed the magnetic field and relativistic 
leptons coexist uniformly throughout the jet. If, alternatively, the magnetic field is enhanced at 
and beyond knot A, or is strong only in localized structures (say a thin surface layer around the 
jet, or the knots within the jet), the inner jet may still have Pb "C Pel throughout most of its 
volume, in accord with our models. 



5.3. Jet Power 

We can also compute the total jet power (energy flux) in terms of the pressure at HST-1. In 
Section 13.11 through Section 13.41 we worked with the total energy flow down the jet (Equations [3] 
or I18p . However, when considering the power the jet can deposit in its surroundings, the "usable" 
energy (kinetic plus internal, omitting rest mass energy of the plasma particle s) is more german e. 
Thus, we write the jet power in the absence of significant magnetic energy as (|Owen et alJl2000l ) 



E, 



(7, ~ 1) 

7i 



Pj c + 



—Pi 
- 1 3 



(22) 



The first term in square brackets is the "useful" kinetic energy of the flow; the second term is the 
internal energy. This equation can also be written as 



(7i 



1) 



+ 1 



(23) 



Because our jet models are mildly relativistic ('jj ~ a few), and have Xj ~ 0(1), both kinetic and 
internal energies contribute significantly to the total jet power. 

Let us evaluate Equation (f23j) at HST-1, where rj ~ 1.3 x 10 19 cm and Pjo = 10~ 8 Po 5 s dyn 



For our fast jet, with jjq = 7.5, Xjo = 0.11 and Tjq = 1.63, we get Ej 



2.0 P .8 x 10 44 erg 

l43 



-i 



cm 

\ For our slow jet, with 7j = 4.4, Xjo = 0.08 and Tj = 1.64, we get Ej ~ 8.2 P ,g x 10 4 " erg 
These values are more or less consistent with estimates for the M 87 jet power in the literature, 
2 x 10 43 erg s" 1 , required to drive the weak shocks associated with the 



which range from Ej 

present jet activity ( Forman et al. 20051 ). to Ej > 10 44 erg s" 1 , based either o n the energy required 
to ex pand the inner lobes into the surrounding ISM over a ~ 1 Myr lifetime (IBicknell &: Begelman 
1993), or on the minimum jet power required at knot A as derived from synchrotron theory (as in 
Section K2 \ lOwen et al.l feOOtf l). 
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6. A Pseudo-Synchrotron Image of Our Model Jet 

In this section we illustrate the effects of the KH modes on the appearance of a jet by producing 
a model image which we compare to the real M87 jet. We choose the fast-intermediate-case jet 
model as being the most likely scenario for the real jet. We compute the intrinsic pressure and 
velocity structure along a jet with a mix of helical, elliptical, triangular and rectangular surface 
modes, chosen so that the elliptical surface mode dominates the structure. We use these results to 
generate a pseudo-synchrotron intensity image that includes all differential Doppler boosting and 
light travel time effects associated with the relativistically moving pressure and velocity structures. 
We view this model jet as being "representative" of a jet like that in M87, but have not attempted 
to create an exact match to the real jet. 

6.1. Pressure and Velocity Structure of the KH Modes 

We compute the pressure and velocity flow fields associated with relativistically moving normal 
mode structures in the limit of weakly non-linear, adiabatic compressions. To do this, we need 
expressions for the spatial structures of the perturbations associated with each KH mode. In the 
linear regime the displacement amplitude of the jet surface associated with a normal mode obeys 



A(z) = Ajo exp 



\ki{z)\dz 



(24) 



where Ajq is the amplitude at Zjo of a sinusoidal displacement to the jet's surface, ki{z) is the 
spatial growth rate for that mode (Figure [7]), and f\kj(z )\dz is " the n umber of growth lengths 
along the jet." We use the expressions found in Section 3 of lHardeel ( 2000i ) to calculate the pressure 
and velocity fluctuations. These grow more slowly than the surface displacement, according to 
A(z)/rj(z) rather than A(z), because the displacement amplitude as a fraction of the jet's radius 
governs the magnitude of velocity and pressure fluctuations. Thus, for constant jet expansion, 



A(z) Aj 



jo 



exp 



Iz jQ h{z)dz 



rj(z) r j0 [l + (z- Zjoty/rjo] 



(25) 



where rj = rjQ + {z — Zjo)ip and tp is the half opening angle of the jet. The small opening angles 

15° or 25° is the viewing angle), are significantly less 
~ 1 at HST-1, as is required for self-consistency. The 



in our models, ip ~ 0.0523 sin 6 (where 9 
than the relativistic Mach angle, ~ (jjMj 
spatial growth length, £ = \ki\~ 1 , for the elliptical surface mode between HST-1 and knot A when 
rescaled to the jet's radius at HST-1 increases from £ ~ lOrjo to £ ~ 20rjQ. For a typical value of 
£ ~ 15rjo ~ 60 pc, pressure and velocity fluctuations grow by an order of magnitude in ~ 150 pc. 

We have explored a number of different mixes for four lower-order KH modes (helical, ellip- 
tical, triangular, rectangular) using spatial growth rates for the fast-intermediate-case jet model 
to evaluate Equation ([25]) . Our goal is not to reproduce the M87 jet exactly but to generate the 
most educational example. To do this we have arbitrarily capped helical and elliptical modes at 
an identical amplitude so that the pressure fluctuation induced by the helical mode is about 50% 
of that induced by the elliptical mode. We also cap the triangular and rectangular modes at a 
smaller amplitude, to create a pressure fluctuation also about 50% of that induced by the ellipti- 
cal mode. With these amplitudes the combined modes produce a maximum pressure fluctuation 
comparable to the pressure fiel d in the unperturb e d jet, i.e., AP S.Pj- This is consistent with 



simulations ( Hardee et al.l fl998l : lAgudo et al.ll200TI : lHardee &: Hugheal2003l ) which show that the 
linearized equations satisfactorily model the pressure and velocity structure in the simulations at 
this fluctuation level. 
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Our results are shown in Figure [8l The top panel indicates how rapidly a small initial pertur- 
bation with amplitude A(z)/rj(z) = 0.01 at HST-1 would grow to our amplitude caps on the M87 
jet for frequency-resonant helical, elliptical, triangular and rectangular surface modes. The lower 
panels of Figure [8] show one-dimensional pressure and velocity cuts, made in a Cartesian coordinate 
system at y = and x/rj = 0.22, 0.44, 0.66, 0.88. Thus, v x = v r is a "radial" velocity component 
in cylindrical coordinates and v y = is a toroidal velocity component. 
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Fig. 8. — Top panel shows the amplitudes, A/rj, for (Hs) helical, (Es) elliptical, (Ts) triangular & (Rs) rectangular 
surface modes used in the lower panels to compute intrinsic pressure and velocity fluctuations. Quantities are 
plotted against z m , the intrinsic distance down the jet. In pressure and velocity panels ID cuts are shown at: 
r/rj = x/r-j = 0.22 (red dashed lines), 0.44 (green dashed lines), 0.66 (blue dashed lines), 0.88 (black solid lines). 

The results illustrate the complex structure that can result from multiple mode superpositions. 
The dominant wavelength of the pressure fluctuation, increasing along the jet but typically ~ 300 
pc, is the wavelength of the elliptical surface mode corresponding to a 180° rotation of the elliptical 
distortion to the jet's cross section. Variation in the fluctuations comes from beating between the 
different wavelengths and wave speeds associated with the other KH modes. The strongest pressure 
fluctuations occur near the jet's surface. The velocity fluctuation in the axial direction, j z v z , is 
small. The modest increase in "radial" velocity, v r , towards the jet surface indicates jet expansion. 
There is more fluctuation in radial v r and toroidal motion near the jet's surface. 

6.2. The Pseudo- Synchrotron Model Image 



In order to image our model jet, w e create a three-di mensional, Cartesian, data cube which 
we use in a modified version of RADIO ( Clarke et al. 19891 ). Each grid element in the data cube is 
derived from the density, pressure and velocity values we calculate throughout the jet, as in Figure 
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[8l These values are suitably corrected for all light-travel-time effects associated with relativistically 
moving structures seen at a viewing angle of 15°. We use a pseudo-synchrotron emissivity, based 
on standard synchrotron theory for a power-law distribution of relativistic lepton energies, with 
Doppler boosting included: 

>2a B l+a S 2+a u -a _ ( 26 ) 



e v oc p)- 2a Pi 



Here, e v is the volume emissivity per Hz; the spectral index a = 0.55 (appropriate for the M87 
jet in the radio band; Frazer Owen, private communication 2009) is taken as constant. This 
expression assumes the leptons are coupled to the total mass density in the je t, are subje c t only 
to adiabatic compression, and ignores radiative losses and in situ energization. iJones et al.l ( 19991 ) 



have shown that Equation (|26|) provides an acceptable image of jet structure in many cases when 

2/3 

the relativistic particles are not tracked explicitly. We assume the field strength obeys B oc p- , for 

a passive, disordered magnetic field. We take Pj(z) oc pj(z) ' 325 , to describe the pressure decline 
in our intermediate-case model jet, and on average the emissivity drops according to 6^(z) oc 
Pj(z) 1 ' 3 oc -Pj(z) 4,0 . Adiabatic pressure fluctuations obey P oc p r around the local values, where 
4/3 < T < 5/3 (equation [6|) . Thus, a factor ~ 1.5 — 2 pressure increase towards the outer edge 
of the jet in a filament (Figure [8|) increases the local emissivity, relative to the mean value, by 
e v oc pj 1 - 1 ™ 3 )^ „ pi-7 _ P i.8 ) thus by a factor „ 2 _ 3 . 

We present our model image in Figured The upper panel shows the basic image, with intensity 
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Fig. 9. — Pseudo-synchrotron intensity image (upper). The axes are in units of arcsec and contours are at the levels 
0.1, 1.25, 2.5, 5, 7, 10, 14, 20, 28, and 40 intensity units per pixel. The resolution is A = r j0 /2 = 0.0262", a pixel is 
a A x A square, and four pixels span the jet diameter at 1". Transverse intensity profiles (lower), taken across the 
jet; axes for these profiles are (vertical), intensity units per pixel and (horizontal) arcsec across the jet. Each profile 
is labelled with its position along the jet in arcsec. 

levels indicated by color scale and contour levels. The lower panels show slices across the model jet 
at ten locations along the jet. At a fixed spatial location, the filaments would rotate temporally 
in the clockwise direction when viewed towards the central engine by an observer, corresponding 
to a clockwise spatial twist of the filaments when viewed from the central engine outwards. If 
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the opposite twist were used the transverse asymmetry seen in Figure [9] would be reversed. The 
importance of relativistic aberration and jet orientation is illustrated by comparing Figure [3] to 
Figure [9l The jet in Figure El chosen to illustrate the intrinsic nature of the KH elliptial mode, 
is at rest in the plane of the sky. The jet in Figure [9] has a similar intrinsic pressure structure, 
dominated by the elliptical mode, but it is moving relativistically and sits at 15° to the line of sight. 
The differences in its appearance, compared to the jet in Figure [3l are due to its speed, orientation 
and the presence of additional dynamically significant KH modes. 

6.3. Comparing the Model Image to the Real Image 

While our model jet image is not tuned to reproduce the M 87 jet, we gain considerable insight 
by comparing our model jet to the 15-GHz image of the M 87 je t in Fi gure [2] of this paper (also 
Figures 5, 6 and 7 of OHC), and the NUV image of iMadrid etaD (120071 ). Our model image shows 



some of the features seen in the real image, but differs revealingly in others. 

Along the model jet the average intensity drops by a factor ~ 10 between 1" and 12". This 
drop is much less than might be expected for an expanding jet which does not decelerate. In a 
constant-speed jet, expansion by a factor ~ 12 between 1" and 12" would cause the density to drop 
by a factor ~ 140, and the pressure to drop adiabatically by ~ 4000, leading to a drastic decline in 
the synchrotron emissivity. Our intensity decline is still somewhat more than in the radio image of 
the real M87 jet; in Figure [2] the jet intensity measured between the bright knots drops by a factor 
~ 2 between HST-1 and knot A. Possible reasons for this difference not accounted for in Equation 
(26) include a slower decline with expansion in an ordered magnetic field and continuous in situ 
energization near to the jet's surface decoupling the radiating leptons from the dynamical density 
and pressure field. These combined effects could lead to the observed slower overall intensity decline 
and to the observed limb brightening which requires an emissivity from the outer 20% of the jet 
~ 5 times higher than that of the inner part of the jet (see OHC). 

Our model image and the accompanying slices convey an overall impression of filaments wrap- 
ping around the jet. The twisted filaments can also be traced in the slices in Figure [9l The four 
slices from 1.7" to 3.75" show winding filaments which brighten alternately at the top and bot- 
tom edge of the jet, but followed by rapid spatial change from 3.75" to 3.85". Various different 
transverse structures are illustrated by the slices from 5.0" to 10.45" including an overall decline in 
intensity and lack of sharp filament structure beyond 8.0". For comparison, filaments in the real 
jet are narrower than those in our image and cross the jet at a more oblique angle, e.g., profiles 
across the real jet (in OHC), from ~ 6" to ~ 10.5". 

The model image is much more asymmetric past ~ 4" than the real jet, with the upper half of 
the model image being ~ 2 — 3 times brighter than the lower half. In large part this is a result of one 
elliptical filament being much brighter than the other. The brighter elliptical filament can easily 
be traced in the intensity image in Figure [9] as it has produced bright spots along the jet top. This 
filament completes a 360° wrap between 1.25 — 3.25", again between 3.25—5.5", yet again between 
5.5—8", and finally again between 8—11". The filament twists around the jet with wavelength varying 
from 2" — 3" along the jet in agreement with the results from LHE. The second elliptical filament 
comes over the jet top about midway between the bright spots associated with the first filament 
and no bright spots are apparent. The difference between the two filaments occurs because the 
pressure structure accompanying the included helical mode adds pressure to one elliptical filament 
and subtracts from the other. In the real jet the two filaments are comparable and no bright spots 
are seen along the top edge. 

The model image leads us to conclude that the helical mode on the real jet must be at very low 
amplitude and dynamically insignificant inside knot A in order for both filaments to be comparable. 
The model image requires a somewhat lower elliptical amplitude in order to reduce the bright spots 
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along the jet top that are not seen on the real jet. In our model image low- level triangular and 
rectangular modes add some complexity, as expected from Figure [HJ and make it more similar to 
the real jet. Note that our model emissivity is assumed to have no radial dependence other than 
that caused by the KH instability. The diminished emissivity in the inner 80% of the real jet would 
also lead to differences between the model and real jet images. 

Even with elimination of the helical mode and reduced amplitude for the elliptical mode, our 
model image would still show an asymmetry. The asymmetry is not due to flow around the helix; 
Figure [8] shows that transverse velocity fluctuations, v r ^^,0.1c, are too small to produce significant 
differential Doppler boosting. Rather, the asymmetry accompanies any helical pattern projected on 
the sky plane. The amount of asymmetry depends on the speed of the helix and its radius relative 
to its wavelength. The helix appears symmetric if viewed close to the critical angle, 8 C = cos -1 /3 W , 
but becomes progressively more asymmetric as the viewing angle moves away from 6 C (Perucho 
et al., in preparation). In the model image the asymmetry is minimal in the inner 4" where the 
filament winds around the jet surface and r,- <C \ ob /2ir. However, with our viewing angle 6 <C C , 
the asymmetry becomes evident at larger distance when, rj < X ob /2tt and cusps appear along the 
jet top beginning at about 8". There is some indication of this behavior on the real jet beyond 
knot F. At about the position of knot A our helical filaments would loop back on themselves at the 
top of the jet and this behavior might be indicated by the complicated structure of knot A in the 
real jet (see Figure 2). 

The real jet contains bright knots inside the jet edges, which are not associated with filament 
crossings, and appear to be features internal to the jet rather than on its surface. The knots in 
the real jet are enhanced in the radio image by a similar factor of ~ 3 — 5 compared to the fainter 
interknot background, but are relatively brighter and more compact in the optical image. Our 
model image shows that these bright knots in the real jet cannot be the result of twisting KH 
modes. The apparent knot location inside the real jet, and the fact that details of knot structure 
change with the observing band, suggest that an additional mechanism is operating in the jet 
interior. 

7. Summary 

We have presented a model in which the dual, helically twisted filaments in the inner M87 
jet are caused by the Kelvin-Helmholtz (KH) instability. We based our model on the qualitative 
similarity between the observed filaments and the likely development - linear growth and nonlinear 
saturation - of the KH elliptical surface mode. The model was developed in three stages. We first 
constructed baseline, steady-state jet models. We then carried out a KH analysis of the baseline 
models. Finally, we created a pseudo-synchrotron image of KH modes in a representative model 
jet, and compared it to the real jet in M87. In this section we summarize the key results of each 
step of our analysis. 

7.1. A Steady-State Jet Model 

Proper motions of bright features in the jet are consistent with a faster jet flow moving through 
a slower wave pattern. Based on these data, we modelled the (unperturbed) inner jet as a weakly 
magnetized, relativistic, steady-state flow, oriented at 15° — 25° to the line of sight, and found the 
following. 

The jet flow decelerates between HST-1 and knot A. It begins with 7^0 ~ 4 — 8 at HST-1, and 
slows to 7j ~ 2 — 3 by knot A. This deceleration is (almost) directly observed in the superluminal 
proper motions of the faster bits of the bright knots. The range of jj values reflects variation in 
reported proper motions and uncertainty in the viewing angle. 
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The wave (pattern) accelerates along the same part of the jet. This is required by the observed 
increase in filament wavelength along the jet (from LHE), and is consistent with subluminal proper 
motions detected in slower bits of the bright knots. 

The jet plasma is internally hot but subrelativistic. If the jet conserves mass and energy as 
it decelerates, its internal energy must be subrelativistic with Pj ~ .04pjC 2 at HST-1 and heating 
to Pj ~ 0.5/jjC 2 by knot A. A modest ( <, 20%) loss of energy from the jet flow to a surrounding 
cocoon does not significantly change this result. Thus, although the jet is detected by means of 
the highly relativistic leptons (P rc \ S> p r cic 2 ) responsible for its synchrotron radiation, its internal 
energy is dominated by a substantially (P t h ^5 PthC 2 ) cooler component. 

7.2. KH Stability Analysis 

We analyzed the KH instability of our steady-state jet models. Our analysis requires that the 
jet be in pressure balance with a surrounding "cocoon" . We determined what conditions must exist 
in this cocoon in order to account for the observed filaments in terms of the KH elliptical surface 
mode without any assumptions about how close to resonance the KH elliptical mode might be. Our 
results from this analysis are as follows. 

The jet is surrounded by a hot cocoon. The cocoon must be denser and cooler than the jet, but 
still quite hot compared to the galactic ISM. At HST-1 the cocoon must have P c ~ 0.002 — 0.003p c c 2 ; 
its mass density must be ~ 10 — 20 times that of the jet. By knot A it must have P c ~ 0.4 — 0.7 p c c 2 , 
and its mass density must be comparable to that of the jet. These results are not derived from 
thermal modeling of the cocoon, but rather are required (through the solution of the KH dispersion 
relation) in order to match the increase in filament wavelength observed along the jet. 

The filaments are KH instabilities close to resonance. The wavelength of the filaments and the 
estimated wave speed are found to be close to the resonant (fastest growing) wavelength expected 
for the KH instability under the conditions we derived for the jet and the cocoon. This was not 
assumed a priori. Because we expect random perturbations to couple to the system's resonant 
frequency, we take this as circumstantial evidence supporting our hypothesis that the filaments in 
the M 87 jet are produced by KH instability. 

The filaments are saturated KH elliptical modes. We determined growth rates for the KH 
modes, and verified that they are rapid enough for the observed filaments to have developed by the 
position of HST-1. Because our model is based on a linear analysis, we cannot analytically predict 
the level at which the KH modes will saturate. However, we suggest the elliptical surface mode is 
responsible for the filaments, because that mode is the most likely to dominate in the nonlinear 
stage without disrupting the jet. The lack of any detectable helical twisting between the core and 
knot A requires that any helical perturbation to the jet be very small. 

7.3. Synchrotron Image of the Model Jet 

We created a pseudo-synchrotron image of one of our fast-jet models, which was modified 
to include helical through rectangular KH surface wave modes. The image took light-travel-time 
and relativistic wave motion into account for our assumed viewing angle. We did not attempt to 
match specific features in the real jet, but comparison between our model image and the real image 
revealed the following. 

The surface brightness along the model jet declines more rapidly than in the real jet but much 
less so than in adiabatic expansion. This is directly due to our assumption that the jet plasma de- 
celerates between HST-1 and knot A and energy flux is conserved. The deceleration and conversion 
of flow to internal energy offsets the radial expansion of the jet, and keeps the psuedo-synchrotron 
emissivity high. The less rapid decline in the real jet indicates that additional in situ processes 
must be operating. 
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The twisted filament structure in the model jet differs from that in the real jet. The filaments 
in the real jet appear more symmetric and are narrower than those in the model jet. On average 
the model jet filaments are brighter along one edge due to the angle at which the helically twisted 
filaments are viewed. One of the two filaments in the model jet has been overly enhanced by the 
inclusion of the helical surface mode. This indicates that any helical perturbation to the jet inside 
knot A must be very small or nonexistent. The real jet is more limb brightened, especially in the 
radio, than is the model jet. This suggests that additional processes are involved in the appearance 
of the real jet on the edge of the flow. 

The bright knots in the real jet are not produced in the model jet. The bright knots in the 
model jet that lie along one edge are due to the angle at which the helically twisted filaments are 
viewed and to the inclusion of the helical surface mode. The bright knots in the real jet tend to 
have a more complex structure than those in the model jet, and also tend to lie inside the jet. This 
indicates that additional processes not associated with KH instability are involved in the real jet 
inside the flow. 



8. Discussion 

How might our model fit within the larger picture of what is known about the inner M 87 jet 
and its environment? Specifically, let us consider the presence of a hot cocoon, the origin of the 
knots, and the disruption of the jet past knot A. 

8.1. Is There an Invisible Cocoon? 

A key part of our analysis is the existence of a hot, tenuous, high-pressure cocoon surrounding 
the inner jet in M87. Such a cocoon must exist if the filaments in the jet are caused by the KH 
instability. Specifically, our models require the cocoon to satisfy P c ~ 0.002 — 0.5p c c 2 , which is 
substantially hotter than the T ~ 10 7 K typical of the inner ISM. Our specific example, in Section 
[5l suggests a cocoon density ~ 10 times lower than the inner ISM at HST-1, and ~ 10 3 times lower 
than the inner ISM at knot A. 

At first glance this requirement seems to contradict the data, because no such cocoon has been 
detected; but it may be hard to see. Any X-ray detection of the cocoon would be indirect; it is 
too tenuous, and possibly too hot, to be X-ray bright. It could appear as a "hole " in an X-ray 
imag e, but the wealth of loops, filaments and cavities in the inner X-ray region ( For man et al.l 



2005) could easily hide such a structure. Direct detection of the cocoon in the radio is equally 



difficult, because uncertainties in relativistic particle energy and magnetic field in the cocoon make 
predictions of its radio brightness elusive. Figure [1] shows that the observed radio halo is relatively 
faint around the inner jet. This could mean that the region immediately around the jet is more 
tenuous than elsewhere in the radio halo, or has a weaker magnetic field, or that the radiating 
leptons have suffered significant losses. These possibilities are consistent with our hypothetical 
cocoon, but by no means confirm its existence. A hint of a low-density cocoon may come from the 
Faraday rot a tion, which is substantially lower along the jet than elsewhere in the inner radio lobes 
Owen et al.l ( 1990I ). While this could be simply due to projection if the jet lies in front of most of 



the inner ISM, it could also indicate a low-density region local to the jet. 

Another possible concern comes from the requirement that the cocoon be at a higher pressure 
than the local ISM. Interpreting the filaments as being KH generated requires the cocoon to be 
in pressure balance with the jet; but the minimum pressure of the jet, P m i n , from synchrotron 
analysis, is higher than the typical pressure of the inner ISM. It follows that the jet-cocoon system 
cannot be in a steady state; it must be dynamic, and expanding outwards into the local ISM. This 
seems consistent with the evidence that the ISM immediately surrounding the jet is disordered and 
turbulent (as in Section [2]). The jet itself is clearly the driver for the inner ISM. It is depositing 



-29 - 



mass and energy directly into the inner radio halo (as in Figure [T]) and our models suggest it can 
be shedding significant energy into its immediate cocoon. 

Yet another concern might be the origin of such a cocoon and why it is so much hotter around 
knot A than at HST-1. While our model does not address this issue, we note that the relativistic 
plasma from the jet is likely to have mixed with the ambient ISM. Such mixing is directly suggested 
by the existence of the radio halo (as in Figure [T]), and is consistent with our example of jet and 
cocoon composition in Section [5j We speculate that relativistic plasma from the jet is mixed 
with the ambient plasma as the jet propagates outward through the galactic ISM. Over time, the 
changing spatial location of the primary mixing and heating region might lead to the required 
temperature and sound speed increase in the cocoon from HST-1 to knot A. 

8.2. What Causes Knots and Limb Brightening? 

There are two important differences between the real jet and our model image. The real jet 
displays the well-known bright knots, which appear to be internal to the jet (as seen in projection) 
and not coincident with filament crossings. Our model jet does not have such complex bright knots. 
In addition, the real jet is limb brightened when seen in the radio, but our model jet is not. 

These differences between the real jet and the model suggest that the jet plasma is being 
energized in situ, in the bright knots as well as at the edge of the jet. This is not a new argument. 
The radiative lifetimes of relativistic leptons emitting X-ray synchrotron emission are too short for 
the particles to survive t he transit from the c ore out to knot A. This result comes directly from the 
strong galactic starlight ( Stawarz et al.|[2005l ). and holds for any plausible magnetic field in the jet. 



It follows that the X-ray-loud leptons must undergo local reacceleration throughout the jet. We 
think it likely that lower-energy particles radiating in optical and radio, are also locally energized 
by the same processes. Spectral behavior along the jet (e.g. Perlman & Wilson 2005, also F. Owen, 
private communication 2009) requires that the acceleration mechanism, whatever it is, maintain an 
approximately constant lepton spectrum throughout the jet. 

One obvious possibility for in situ energization is the KH instability itself. We showed in 
Section 14.31 that the wave speed is supersonic throughout the jet. Thus, the twisted filaments 
provide a continuous pair of twisted shocks, and possible particle acceleration near the jet surface, 
from HST-1 to knot I. In addition, velocity shear at the edge of the jet may amplify the magnetic 
fields there, locally enhancing the synchrotron emissivity close to the edge. If this is the case, the 
weaker limb brightening in the optical, as compared to the radio, must be due to a relatively low, 
high-energy cutoff of the particle distribution near the jet surface. 

However, edge effects are not likely to be the full answer. In general, the knots appear located 
within the jet, not close to its edge (although the HST-1 flaring region was at the jet edge), and 
contain compact, bright emission regions when seen in the optical or NUV. Thus, the data strongly 
suggest that energization internal to the jet is also an important part of the picture. This is 
unsurprising, because the jet is supersonic from HST-1 to knot I, and very likely to contain internal 
shocks. Shocks c ould arise because t he jet flow is unsteady; this is evident from the variability 
of the knots, e.g., iBiretta et al.l ( 19991 ). Shocks could also arise as the jet propagates through the 



turbulent ISM in the galactic core. Particle acceleration at internal shocks could very plausibly 
create the bright knots seen in the radio, optical and X-ray. 

8.3. What Happens Past Knot A? 

We can use our model to speculate on the behavior of the jet downstream of knot A. Radio 
images (see Figure [1]) show the jet narrows, and its centerline oscillates, between knots A and C. 
Our models find the jet is hot and relatively slow at knot A, with sound speed aj ~ 0.5c and Lorentz 
factor 7j > 2.5. If knot A is a shock, the sound speed would be higher and the Lorentz factor would 
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be lower immediately post-shock. In a linear analysis, at a fixed pressure fluctuation, the transverse 
fluid displacement associated with a KH mode is larger on a w eakly superso nic, weakly relativistic 
jet than it is on a highly supersonic and highly relativistic jet (|Hardeell2000l ). Thus, the narrowing 
of the jet past knot A where the flow could be weakly supersonic and weakly relativistic might be 
associated with significant elliptical distortion of the jet's cross section associated with the elliptical 
surface mode of the KH instability. 

In addition, the highly collimated jet flow disrupts significantly downstream of knot C. As we 
noted in Section \2.2\ the KH helical mode does not saturate but grows in amplitude until the jet 
decollimates. We speculate that a slowly growing, long-wavelength helical mode is what disrupts 
the jet downstream of knot A. To test this idea, even with a discontinuity at knot A, we can 
estimate the downstream development of the KH instability in terms of flow conditions at knot A. 
On a weakly supersonic jet, KH stability analysis shows that the supersonic resonance moves to 
very high frequencies and very short wavelengths. In this situation the development of KH surface 
modes can be estimated from the solution of Equation (|A5p . using lujj = w#/2. Our models at knot 
A (e.g., Table 2) have ojrj/uj ~ 1, also pj/p c ~ Xj/Xc ~ 1, 7c ~ 1 and jj£2. With these numbers, 
the growth length for the helical surface mode, £(ujh) = h)\~ 1 ■, is no more than a few jet radii. 
Noting that r, ~ 50 pc at knot A, a helical perturbation there will grow by an order of magnitude 
in ~ 200—300 pc. For comparison, the distance from knot A to knot C is ~ 1.2 kpc with the jet at a 
15° viewing angle. Thus, our results suggest that destabilization and decollimation of the M87 jet 
beyond knot A is determined by conditions in the jet before it reaches that point. With minimal 
energy flux lost from the jet, the decelerating jet flow must heat up; the cocoon must also heat in 
order to account for the observed helical filaments. These conditions lead to rapid growth of the 
KH helical mode, and consequent jet disruption, past knot A. 

Our goal in this paper has been to uncover properties of the M87 jet and its surroundings by 
assuming the KH instability is responsible for the twisted filaments seen in the j et. It is illuminatin g 
to compare our conclusions here to those of a similar analysis of the 3C 120 jet ( Hardee et al. 20051 ) . 



Bright knots and inter-knot emission in the 3C 120 jet on scales of hundreds of pc, suggest conical jet 
flow and a more slowly moving short wavelength helically twisted filament which is consistent with 
a short wavelength helical surface mode. The observed behavior of this mode in 3C 120 requires a 



coolin g, modestly accelerating flow along the expanding jet. Under these conditions, lHardee et al 



( 20051 ) argued that the 3C 120 jet is stabilized to KH disruption by saturation and the declining 



growth rate of the short-wavelength mode along the jet. Their result explains how the 3C 120 jet 
can remain highly collimated to hundreds of kpc from the central engine. By contrast, in this paper 
we find a different result for the M87 jet. We suggest that heating of the M87 jet as it decelerates 
allows the KH helical mode to grow significantly downstream of knot A, leading to jet disruption 
only a few kpc from the central engine. 
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A. Propagation and Growth of Unstable KH Modes 

The propagation and growth of the normal modes of the KH instability depend on the jet 
and cocoon speeds, the thermal state of the jet and the cocoon, the rate of jet expansion, and the 
frequency of the wave relative to a resonant frequency u* and wavelength A* = 2ttv* u] /uj* at which 
the growth rate is a maximum. In this Appendix we review the basic structure of our KH analysis 
and a nalytic approx imations for mode growth below, at, and above the resonant frequency. We 
follow lHarded (120071 ) and refer the reader to that reference for more detail. 

In cylindrical geometry (with jet flow along the z-axis) a random perturbation to the velocity, 
pressure or density can be described in terms of Fourier components of the form 

f(r, (j), Z, t) = f(j.y[Kkz±n<t>-u>t)\ _ ( Al ) 



This preprint was prepared with the AAS IATjrjX macros v5.2. 



-33- 



Here, n is a dimensionless azimuthal wavenumber; n = 0,1,2,3,4... describes pinching, helical, 
elliptical, triangular, rectangular, etc. distortions to the jet cross section. Formall y there are a n 
infinite number of solutions for each "normal mode" n with different radial structure (|Hardeell2000l ) . 
In this paper we only consider the low-n modes, and argue that the elliptical mode dominates the 
structures seen in the M 87 jet. 

To derive the dispersion relation, we assume the jet has a "top hat" profile. We take density, 
pressure and velocity as uniform across the jet, and assume that the external cocoon is in pressure 
balance with the jet, but can otherwise have density and velocity different from those in the jet. 
In this s ituati on, the normal modes are governed by the wave dispersion relation, Equation (2) in 
Harded l|2007l ) : 

CjJUQrj) _CcH^'(C c r j ) (A2) 
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Here, J n and Hn are Bessel and Hankel functions; the primes denote derivatives with respect 
to their arguments. The subscripts j and c refer to the jet and cocoon, and rj is the jet radius. 
In the unmagnetized limit, the terms in Equation (|A2j) are as follows. The term £ is defined as 
£ = "f 2 Wvj 2 , with the plasma energy content W = p(l + x) defined in Equation Q and the specific 
enthalpy \ defined in Equation ([5]). The factor w is defined by w 2 = (u — ku) , for a flow speed 
u = j3c with a Lorentz factor 7. The term Q is defined by £ 2 = j 2 a 2 [w 2 — K 2 a 2 ) where a is the 

sound speed, given by Equation f)15[) . and k 2 = (k — wu/c 2 ) 2 . Formally, this dispersion relation is 
for a cylinder. iHarded (| 19821 ) showed that a similar dispersion relation could be derived for conical 
expansion at much less than the Mach angle. The local dispersion relation above delivers the same 
result if incremented along a conical jet using the approp riate instanta neous local conditions and 
spatial growth becomes power-law instead of exponential ( Hardee 19871 ). 
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For a given jet and cocoon flow speed, solutions of Equation (|A2p depend only on the thermal 
state of the jet and cocoon, which are described by Xj an d Xc respectively. At first glance, Equation 
(|A2|) appears also to depend on the jet and cocoon densities, pj and p c , through the Wj and W c 
terms in the denominators of each side. However, because we assume pressure balance between the 
jet and the cocoon, the ratio Wj/W c in Equation (lA2j) depends only on the ratio of the jet and 
cocoon enthalpies, Xj/Xc- 

The solution of Equation (|A2p relates the wave frequency, u), to a complex wavenumber, k(cj) = 
k r (u) + iki(ui), for a given set of jet and cocoon parameters. Formally the solution to the dispersion 
relation provides an inverse complex phase velocity as a function of oj. Thus, the real part of uj/k 
determines both the wave speed, v w (uS) = [u/k(ui)] r , and the wavelength, X w (oj) = 2itv w (uj)/uj. 
Spatial growth or damping of the wave is determined by the imaginary part of the inverse, i.e. by 
ki/u. A negative value for ki(u) indicates growth with e- folding length, 1(uj) = |/cj(cj)| _1 . 

In general, Equation (lA2j) must be solved numerically, as we do in this paper. However, the 
important propagation properties of this mode can be revealed by analytical approximations. The 
maxim um growth rat e for the elliptical surface mode occurs at a resonant frequency (from Equation 
(16) in IHarded (|2007l )) 
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where the angular frequency u refers to a perturbation propagating along the jet's surface at fixed 
cylindrical coordinate <fi. For the elliptical mode this frequency and the accompanying wavelength 
correspond to a 180° rotation of the elliptical cross section and high-pressure filaments and not 
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a 360° rotation of an individu al high-pressure filament. The wave speed at resonance is (from 
Equation (15) in lHarded (|2007l )) 
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where j a = (1 — a 2 /c 2 ) -1 / 2 is the sonic Lorentz factor. The denominator on the right-hand side 
of Equation (|A3|) shows a rapid increase in the resonant frequency as the jet becomes transonic. 
This means that the jet speed becomes comparable to the sound speed and a c /v^ increases with 
given by Equation (jA4j). Note that the resonant frequency is inversely proportional to the jet 
radius and directly proportional to the cocoon sound speed provided the cocoon sound and flow 
speeds are small. 

In the low- frequency limit (to <C uj* ) the elliptical surface mode propagates and grows according 
to (from Equation (8) in lHarded (|2007l ^ 



k 1 + r\ 1 + r\ 

where r\ = ^Wjjrf.W c . We see that v w = \uj/k(uj)] r can be large or small, according to whether rj 
is large or small. In general, the resonant wave speed is faster than the low-frequency wave speed 
when 7] is small and slower than the low-frequency wave speed when r\ is large. Cocoon flow, u c , 
can significantly affect the wave propagation speed and wavelength if 77 < 1. From the solution to 
Equation (|A5p we also estimate the growth length, £(uj) = |/cj(a;)| -1 . When Uj 3> u c , which is the 
case in this paper, £(uj)/rj ~ rf^Uj /oorj for this low-frequency limit. 
In the high-frequency limit (u ^> uj*) 



k 1 — cijiij/c 2 
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for the unstable branch of the elliptical mode. Here, a minimum possible wave speed occurs for 
aj = c/\/3; for a relativistic jet with Uj ~ c the intrinsic high frequency wave speed is v w > 0.63c. 
Since the condition i) > 1, necessary for high wave speed in the low-frequency limit, typically 
requires aj ^ c/\/3, we see that the high-frequency wave speed exceeds the low-frequency wave 
speed even in this situation. In general, the wave speed at high frequencies exceeds the wave speed 
at frequencies below resonance. In the high-frequency limit cocoon flow plays no role in the wave 
propagation speed. 



